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Abstract 

Various solutions are displayed and analyzed (both analytically and numerically) of 
a recently-introduced many-body problem in the plane which includes both integrable 
and nonintegrable cases (depending on the values of the coupling constants); in par- 
ticular the origin of certain periodic behaviors is explained. The light thereby shone 
on the connection among integrability and analyticity in (complex) time, as well as 
on the emergence of a chaotic behavior (in the guise of a sensitive dependance on the 
initial data) not associated with any local exponential divergence of trajectories in 
phase space, might illuminate interesting phenomena of more general validity than for 
the particular model considered herein. 

1 Introduction 

a 

Recently it was pointed out [1, 2] that the iV-body problem in the plane characterized by 
the following Newtonian equations of motion possesses lots of completely periodic solutions: 

a: 

N 

f n =ukAf n + 2 ^ (r nm )~ 2 (a nm + a' nm kA) 

m=l; mj^n 



X 



^"nv m ' Trim) ^mvn ' fnm) Tnmvn ' frn) ■ (1.1a) 



Here and below the subscripted indices run from 1 to JV (unless otherwise indicated), the 
number of moving particles N is a positive integer, the iV two-vectors r n = r n (t) identify 
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the positions of the moving point-particles in a plane which for notational convenience we 
imagine immersed in ordinary three-dimensional space, so that f n = (x n ,y n ,0); k is the 
unit three-vector orthogonal to that plane, k = (0,0, 1), so that k A f. n = (— y n ,x n ,0); 

Turn = T<n, ^mi ^nm — 2?* n • T m , (1.1b) 

superimposed dots denote of course time derivatives; uj is a real - indeed, without loss of 
generality, positive - constant, lu > 0, which sets the time scale, and to which we associate 
the period 

T = 2it/uj (1.2) 

(when discussing below specific examples we will conveniently set uj = 2ir, so that T = 1; 
but for the moment let us not do that); and a nm and ot' nm are real "coupling constants". 
Note that this iV-body model in the plane is invariant under both translations and rota- 
tions; as we shall see below, it is Hamiltonian; it features one-body and two-body velocity- 
dependent forces; when the latter are missing - namely, when all the two-body coupling 
constants a nm , a' nm vanish, 

Otnm = ol nm = 0, (1.3a) 

it represents the physical situation of N electrically charged particles moving in a plane 
under the influence of a constant magnetic field orthogonal to that plane ("cyclotron"). 
In such a case of course every particle moves uniformly, with the same period T, see (jl.2j) . 
on a circular trajectory the center c n and radius p n of which are determined by the initial 
data (position and velocity) of each particle: 

r n (t) = c n + p n sin (ut) - k A p n cos (ut) , (1.3b) 

Cn = r n (0) + k A p n , p n = f n {0)/uJ. (1.3c) 

Let us begin by tersely reviewing previous findings [1, 2]. First of all we note that it is 
convenient to identify the physical plane in which the motion takes place with the complex 
plane, via the relation 

r n = (x n ,y n ,0) <^ z n = x n + iy n , (1.4) 

whereby the real Newtonian equations of motion in the plane (1.1) become the following 
equations determining the motion of the N points z n = z n (t) in the complex z-plane: 

N 

Z n — iuiZn ~\~ 2 ^ O-nm^n^m / (%n ^m) : (1.5a) 

m=l; rra^n 

with 

Q-nm — ot nm -\- ia nm . (1.5b) 

Hereafter we always use this avatar, (1.5), of the equations of motion (1.1), and we more- 
over exploit the following key observation [3, 4, 1, 2]: via the change of (independent) 
variable 

z(t) = C(r), r = [exp (tut) - 1] / (iw) , (1.6) 
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the system (1.5) becomes 



N 



C = 2 £ a 




(1.7) 



m=l; m^n 



Here and below the underlined notation indicates an iV-vector, say z = (zj, . . . , zjy), 
£ = (Ci> ■ ■ ■ > Civ) and so on, and the primes denote of course differentiations with respect 
to the independent variable r. Note that the constant u has completely disappeared 
from ()1.7)1 ; nor does it feature in the relations among the initial data for (1.5) and (|1.7|) . 
which read simply 



This of course entails that, to obtain the solution z(t) of (1.5) corresponding to a given 
set of initial data z (0), z (0), one can solve 1)1.7)) with the same set of initial data, see 1)1.8)) . 
thereby determine C( r )) an d then use 1)1.6)) to obtain z{t) (hence as well, via (|1.4|) . the 
solution of the initial-value problem for (1.1)). This possibility - to infer the behavior of 
the evolution in the real time variable t of the solutions z(t) of (1.5) (namely as well of 
the solutions of the physical many-body problem in the plane (1.1) ) from the properties 
of the solutions £ (r) of Q1.7|) as functions of the complex variable r - is, as we will see, 
the main tool of our investigation. Indeed, when the real time variable t evolves over 
one period T - say, from t = to t = T - the complex time-like variable r goes from 
r = back to r = by traveling counter-clock-wise - see 1)1. Bj) - on the circular contour C 
centered in the complex r-plane at i/u and having radius 1/uj. Hence whenever all the 
functions Q n (t) - obtained as solutions of ()1.7|) - are holomorphic, as functions of the 
complex variable r, inside that circular contour C - or, more precisely, in the (closed) 
circular disk C defined by that contour C - the corresponding functions z n (t) - namely 
the corresponding solutions of (1.5) - are nonsingular and completely periodic in the real 
time variable t, with period T, see 1)1.2)1 . (This also entails that, if the only singularities 
of a solution £ (r) of 1)1.7)1 inside the disk C are a finite number of rational branch points, 
then the corresponding solution z (t) of (1.5) is again completely periodic in the real time t, 
albeit with a larger period which is then an entire multiple of T - we shall discuss in detail 
this important point below). 

Now we note that all solutions £ (r) of 1)1. 7)1 . corresponding to arbitrary initial data 
C (0), £' (0) (with the only restriction that these data be nonsingular, namely \C' n (0)| < oo 
and Cn (0) 7^ Cm(0), see the right-hand side of 1)1.7)1 ) . yield solutions C ( r ) which are 
holomorphic in the neighborhood of r = 0, as implied by the standard theorem which 
guarantees the existence, uniqueness and analyticity of the solutions of analytic ODEs, in 
a sufficiently small circular disk D centered in the complex r-plane at the origin, r = 0, 
where the initial data defining the solution are given. The size of this disk D is determined 
by the location of the singularity of ( (r) closest to the origin, and the structure of the 
right-hand side of 1)1.7)1 clearly entails that a lower estimate of this distance - namely of 
the radius p of D - reads as follows: 



with R a positive constant (dependent on the values of the coupling constants a nm but 
not on the initial data) and £ respectively £' providing lower respectively upper estimates 



!(0) = C(0) 



i(o) = c'(o). 



(1.8) 



p > RC/C 



(1.9a) 
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of the moduli of Q n (0) — Q m (0) respectively Q' n (0), say 

C= min |C«(0)-Cm(0)|, (1.9b) 

n,m=l,...,N ; n=£m 

(> = n m^ N \(' n (0)\ (1.9c) 

(for a derivation of this formula, including an explicit expression for R, see Appendix A). 
Let us now assume the initial data, see ljl.8|) . to entail (via (1.9)) that p > 2 /to. Then 
the disk D includes the disk C, and this entails that the solutions Q (r) of Q1.7|) are 
holomorphic functions of the complex variable r in the (closed) disk C, hence (see Q1.6)) ) 
the corresponding solutions z(t) of (1.5) are completely periodic functions of the real 
variable t, with period T, see (|1.2j) . 

This observation, together with the lower estimate (1.9) of p, imply the existence of a 
set, of nonvanishing (in fact, infinite) measure in phase space, of initial data z (0), i (0) 
which yield completely periodic solutions of (1.5) (hence as well of (1.1)). But before 
formulating this finding [1, 2] in the guise of the following Proposition let us interject 
the following obvious 

Remark 1.1. If z (t) is the solution of (1.5) corresponding to initial data z (0), z (0), then 
z (t) = cz (bt) is the solution of the equations of motion that obtain from (1.5) by replacing 
in it uj with uj = buj, and of course the corresponding initial data read 5(0) = cz(0), 
z(0) = bcz(0). Here b, c are of course arbitrary (nonvanishing!) rescaling constants. 

Let us now formulate Proposition !!. 21 (correcting thereby the partially incorrect formu- 
lation of this result given in Refs. [1, 2]). 

Proposition 1.2. Let z{t) be the solution of (1.5) with 

uj = bu, (1.10a) 

corresponding to the assigned initial data 

z(0) = cu, 1(0) = pv, [with u n ^ u m if n ^ m] , (1.10b) 

where the positive numbers b, c, p play the role of scaling parameters (as we shall im- 
mediately see). Then the solution z(t) is completely periodic with period T, see \1. 6 J\) . 



z(t + T) = z(t), (1.11) 

provided one of the following conditions hold: 

(i) for given a nm , uj, z (0) andv, the scaling number p, hence as well the initial velocities 
z n (0), are sufficiently small: < p < p c , where p c is a positive number, p c > 0, whose 
value depends on the given quantities; 

(ii) for given a nm , to, i(0) and u, the scaling number c is sufficiently large, c > c c 
(hence the initial positions of the N particles in the plane are sufficiently well separated), 
with c c a positive number, c c > 0, whose value depends on the given quantities; 

(Hi) for given a nm , uj, z(0) and i(0) ; the scaling number b, hence as well the circular 
frequency uj = buj, is sufficiently large, b > b c , where b c is a positive number, b c > 0, whose 
value depends on the given quantities. 
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Remark 1.3. The first two formulations (items (i) and (iij) of Proposition lTT^l refer to the 
same equations of motion, with modified (rescaled) initial conditions; the third formulation 
(item (in)) refers to different equations of motion (due to the change via rescaling of 
the constant oS) with the same initial conditions. But in fact these 3 formulations are 
completely equivalent, see Bemark ll.ll 

We have thereby reviewed (and rectified!) the results of References [1, 2], to an extent 
sufficient to make this paper self-contained. But before focusing on the specific, novel, 
contributions of this paper three additional facts are now recalled. The first two are 
relevant to the case - to which attention will be hereafter restricted - in which the (possibly 
complex) coupling constants a nm depend symmetrically on their two indices: 

Q"nm — O mn . (1.12) 

Then the equations of motion (1.5) are Hamiltonian, being derivable in the standard 
manner from the Hamiltonian 

N ( N "j 

H{z,p)=J2\ iuJZ n/ c + ex P (cPn) [I ^ ~ Z ^~ a " m > (1-13) 

n=l I m=l;m^n I 

(and let us recall that this entails that the "physical" equations of motion (1.1) are as 
well Hamiltonian [2]). Note the presence, in this expression of the Hamiltonian function 
H(z,p), of the arbitrary (nonvanishing) constant c, which does not feature in the equations 
of motion (1.5). Also note that H(z,p), in contrast to the equations of motion (1.5), is 
not quite invariant under translations (z n — > z n = z n + Zq), although the only effect of 
such a translation on H(z,p) is addition of a constant. 

Moreover, as clearly implied by the equations of motion (1.5) with (|1.12|) . the center of 
mass, 

N 

Z(t) = N- 1 Y,Zn(t), (1.14) 

n=l 

moves periodically (with period T, see p.2|l ) on a circular trajectory (in the complex 
z-plane): 

Z(t) = Z (0) + Z (0) [exp (iut) - 1] / (iu) . (1.15) 

The third and last fact we like to recall is that, if all the coupling constants in (1.5) are 
unity, 

a nm = 1, (1-16) 

then the equations of motion (1.5) are integrable indeed solvable [5, 2], the solution of the 
corresponding initial-value problem being given by the following neat prescription: the 
N coordinates z n (t) which constitute the solution z(t) of the equations of motion (1.5) 
corresponding to the initial data 2(0), i(0) are the N roots of the following algebraic 
equation in the variable z: 

N 

Zm (0) /[z-z m (0)] = iu/ [exp {iut) - 1] . (1.17) 

m=l 
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Note that, after elimination of all denominators, this is a polynomial equation of degree iV 
for the variable z, with all coefficients of the polynomial periodic in t with period T, 
see (|1,2|) . Hence the set z(t) of its N zeros is as well periodic with period T. It is 
therefore clear that, in this special solvable case, see ()1.16(l . all nonsingular solutions 
of (1.5) are completely periodic, with a period that is either T or an integer multiple of T 
(the latter possibility arises because the zeros of the polynomial - which is itself periodic 
with period T - may get reshuffled through the motion - but of course that integer multiple 
cannot exceed N\ - indeed it cannot exceed the combinatorial factor p (N) defined as the 
smallest integer p (N) such that the iteration p (N) times of any permutation of N objects 
yield unity, and of course p (N) <C Nl for large N). 

As we just saw, Proposition 11.21 entails that the equations of motion (1.5) (as well as 
the equivalent "physical" equations of motion (1.1)) possess a lot of nonsingular and com- 
pletely periodic solutions. This finding [1, 2] was (of course!) validated by simulations 
performed via a computer program created by one of us (MS) to solve numerically the 
equations of motion (1.5). The new findings discussed in this paper emerged from the 
(rather successful!) effort to understand certain remarkable features of these numerical 
simulations - in particular the existence, in the cases characterized by real and rational 
values of the coupling constants a nm in (1.5), of additional sets of initial data, also of non- 
vanishing measure in phase space, yielding nonsingular and completely periodic solutions 
with periods which are integer multiples of T, see (jl.2|) - via a mechanism, at which we 
already hinted above, associated with rational branch points of the solutions ( (r), which 
is indeed also responsible in the solvable case (with a nm = 1) to yield completely periodic 
solutions with such larger periods (integer multiples of T). These new findings also suggest 
the existence of other integrable - but presumably not solvable - cases of the many-body 
problem (1.1), as well as of nonintegrable cases for which however integrable behaviors, 
such as those associated with completely periodic motions, emerge from sectors of initial 
data having nonvanishing measure in phase space; and they display as well a mechanism 
for the emergence of a chaotic behavior (for other regions of phase space, of course only 
in the nonintegrable cases) characterized by a sensitive dependence on the initial data 
(the signature of chaos!) which is not associated with any local exponential divergence 
of trajectories over time - in analogy to the type of chaotic behavior that ensues when a 
particle moves freely (except for the reflections on the borders) inside, say, a triangular 
plane billiard whose angles are irrational fractions of ir. 

These findings shall be reviewed at the end of this paper, after they have been discussed 
in the body of it. In particular, in the following Section 2 we discuss the conservation laws 
associated with the equations of motion (J1.7|) ; in Section 3 we exhibit certain similarity 
solutions of these equations of motion; in Section 4 we discuss the solution of the two-body 
problem, in somewhat more detail than it was done in [1]; in Section 5 we investigate 
the analytic structure of the solutions of l|1.7|) . in particular the nature of the branch 
points, in the complex r-plane, featured by the solutions of these equations of motion, 
and the structure of the associated Riemann surfaces, which is crucial to determine the 
behavior - be it completely periodic or chaotic - of the solutions of the "physical" equations 
of motion (1.5) namely (1.1) as the real time t unfolds. The relevance of these findings for 
the motions of the "physical" problem (1.5) or (1.1) are then discussed, and also illustrated 
via the display of numerical simulations, in Section 6. But - to avoid this paper becoming 
excessively long - the treatments of Sections 5 and 6 are restricted to the case in which 
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the real parts of the coupling constants a nm are nonnegative, 

Re(a nm )>0, (1.18) 

a condition which is sufficient (albeit not necessary; see below, in particular (5.8)) to 
exclude the phenomenon of "escape to infinity" ; the general case shall be treated in a sub- 
sequent paper. A final Section 7 outlines how the study of this model can and should be 
pursued. To avoid interrupting the flow of the discourse certain technical developments 
are relegated to four Appendices. 

2 Conserved quantities 

The equations of motion (1.7) are of course (see (1.13)) implied by the Hamiltonian 

N ( N 

= e { ex p n K« - ^' anm \ • 

n=l I m=l;m^n 

Indeed this yields the Hamiltonian equations of motion 

N 

C = dH/dir n = cexp{cir n ) J] [C„ - Cm]~ anm , (2.2a) 

m=l; m^n 

N 

7T' n = -dH/d( n = C- 1 Y, (C + Cm) (Cn" Cm), (2.2b) 
m=l; mj^n 

where to write in more convenient form the second set of equations, (|2.2b|) . we used the 
first set, ()2.2a|) (and we of course assumed validity of the symmetry property (1.12)). 
Clearly logarithmic differentiation of (|2.2aj) yields 

N 

Cn/Cn = c7T n ~ ^ ] a nm (Cn — Cm) / (Cn — Cm) ; (2-3) 
m=l; mj^n 

and these equations, via (|2.2b|) . clearly yields the equations of motion (1.7). 
It is thus seen, from (|2.1|) and (|2.2a|) . that 

N 

H = c~ 1 Y,C (2-4) 

77=1 

hence in this case the fact that the Hamiltonian is a constant of the motion coincides with 
the time-independence of the "center-of-mass" velocity, 

N 

(d/dr)J2C = Q, (2-5) 

77 = 1 

which is itself an immediate consequence of the equations of motion (1.7) (indeed, their 
sum over the index n from 1 to N implies the vanishing of the right-hand side due to 
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the antisymmetry under exchange of the two dummy indices n, m of the summand in the 
double sum - since we always assume (1.12) to hold). 

A second conserved quantity, associated with the translation-invariant character of the 
Hamiltonian 112. IB . is the total momentum II, 



N 



n = ^vr n , 



(2.6) 



n=l 



as it is clear by summing over the index n from 1 to N the equations of motion l)2.2b|) . 

By exponentiation and via H2.2aj) this entails the following version of the second con- 
served quantity: 



N 



n=l 



N 



c n (Cn-ur 

m=l; m^n 



(2.7) 



The fact that this quantity, K, is a constant of motion can also be verified as follows. 
Logarithmic differentiation of (|2.7[) yields 



N 



K'/K = £ 

n=l 

N 

K'/K = Y, 



n=l 



N 

Cn/Cn + a nm (Cn ~ Cm) I (Cn ~~ Cm) 

m=l; m^n 

N 

Cn/Cn ~ 2 ^ 1 a nmCm/ (Cn — Cm 
m=l; m^n 



(2.8a) 
(2.8b) 



The second of these equations, (|2.8b|) . follows from the first, H2.8a() . by exchanging the role 
of the two dummy indices n and m in the first term of the double sum (and of course by 
taking advantage of (1.12)). It is then obvious, see (1.7), that K' vanishes, namely that K 
is a constant of motion. 



3 Similarity solutions 

In this section we report certain special solutions of the evolution equations (1.7), hence as 
well, via (1.6), of the Newtonian equations of motion (1.5) and (1.1). They read (obviously, 
up to arbitrary permutations) 

Cn (t) = c + c n (r - r b f , n = l,2, ...,M, (3.1a) 
Cn(r) = C n , n = M + l,M + 2,...,N, (3.1b) 

with M a positive integer, 2 < M < N, the complex constant arbitrary (it identifies the 
value of r at which this solution has, generally, a branch point; see (j3.1a|) and below), the 
constants C n arbitrary (and largely irrelevant: particles that do not move, see ()3.1b|l . can 
be altogether ignored, since the structure of the right-hand side of (1.7) entails that they 
neither feel nor cause any force), the constant cq also arbitrary (reflecting the translation- 
invariant character of the model) , while the constants T and c n are determined (the latter 
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up to an arbitrary common multiplicative factor) by the following equations, which clearly 
obtain by inserting (3.1) in (1.7): 

r (r - 1) = -r 2 A (3.2a) 



M 

A = -2 a nm c m /(c n - c m ), n=l,2,...,M. (3.3a) 

m=l; mj^n 

The first of these equations, (|3.2ajl . entails 

r = {i + A)~ 1 (3.2b) 

(provided A 7^ —1; if instead A = — 1 the function (r — Tb) r in the right-hand side of IjM.laj) 
must be replaced by exp (/3r), with (3 an arbitrary nonvanishing constant; see below). 

The M equations (|3.3a|) determine A (note the assumed independence of this quantity 
from the index n) and, up to a common rescaling factor, the M constants c m with m = 
1,2, ...,M. The expressions of these M constants cannot be generally obtained (for 
arbitrary a nm ) in explicit form from ()3.3a|) . but this is instead possible for A, hence, 



via (|3.2b|) . for T (provided, as we always assume, the symmetry property (1.12) of the 
coupling constants a nm hold). Indeed by replacing c m in the numerator in the right-hand 
side of (|3.3aj) with (c m — c n ) + c n we get 

M M 

A = 2 ^ Q"nm 2 ^ CL nm C n /{c n c m ), (3.3b) 
m=l;mj^n m=l;m^n 

hence, by summing over n from 1 to M, 

M 

MA = 2 a nm- MA. (3.3c) 

m,n=l; m^n 

To obtain this equation, (|3.3c|) . we exchanged the order of the second double sum over n 
and m in the right-hand side, we took advantage of the symmetry property (1.12), and 
we used again (|3.3aj) . From (|3.3c|) (and taking again advantage of the symmetry property 
(1.12) so as to avoid any double counting in the sum) we finally get 

M 

A=(2/M) a nm- (3.3d) 

n,m=l; n>m 

In conclusion we see that the system (1.7), for any arbitrary choice of the coupling 
constants a nm , possesses the (exact) solution (3.1), with an arbitrary choice of the positive 
integer M in the range 2 < M < N and of the N — M (complex) constants C m (for m 
in the range M + 1 < m < N), and with the M constants c n (for n in the range 1 < 
n < M) determined, up to a common rescaling factor, by the system of M algebraic 
equations ()3.3a|) . where the constant A is given (provided, as we always assume, there holds 
the symmetry property (1.12)) by the simple formula (|3.3d|) . The (generally complex) 
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constant tj in (|3.1aj) is also an arbitrary constant, and its significance is that, at r = r&, 
the M particles labeled from 1 to M all "collide" if Re (r) > or "escape to infinity" if 
Re (r) < 0, yielding of course a singularity of the equations of motion, as evidenced by 
the right-hand side of (1.7). Note that the solutions Cn (t) possess, at this point r = t&, 
a branch point singularity, whose nature is characterized by the exponent T, see (|3.1aj) and 
()3.2b|) with (|3.3d|) . Of course analogous solutions exist for any arbitrary permutation of 
the N indices n. (Note that above we put inverted commas around the word "collide" or 
"escape to infinity" - to underline that this is not, generally, a true "physical" phenomenon, 
since it occurs at some complex value of the r variable; see below). 

The corresponding solutions z{t) of the equations of motion (1.5), hence as well (via 
(1.4)) of the Newtonian equations of motion (1.1), obtain from these solutions C ( r ) 
via (1.6). Therefore when the complex constant r& falls outside the disk C (centered, 
in the complex r-plane, at r = i/u, and with radius l/u>), these solutions z_{t) are non- 
singular and completely periodic with period T, see (1.2), while if the complex point r& 
falls inside the disk C, the solutions z (t) are again nonsingular but generally not periodic, 
unless the exponent T, see (|3.1a|) and ()3.2b|) with l)3.3dl) is rational, V = p/q (with p and q 
two coprime integers, and, say, q positive, q > 1), in which case the solutions z (t) are again 
completely periodic, but with period T = qT (as clearly entailed by ()3.1a|) with (1.6)). In 
the special case in which the point falls exactly on the circle C that constitutes the 
boundary of the disk C, so that the formula t& = [exp {iujt c ) — 1] / (iu>) defines mod (T) 
a real time t c , then the solution z (t) corresponds to a special Newtonian motion in which, 
at the "collision time" t c , M particles simultaneously collide (if Re (r) > 0) or escape to 
infinity (if Re (r) < 0). 

This analysis applies provided A ^ — 1 (see ()3.3d|0 . so that T is well defined (see 
(|3.2b|l ). If instead A = —1, the similarity solutions ( (r) are entire functions of r, since 
the expression (|3.1a() is then replaced by the formula 



where (5 is an arbitrary constant (just as t\, is an arbitrary constant in (|3.1a|l ). Hence in 
this case the corresponding similarity solutions z_ (t) are always nonsingular and completely 
periodic with period T, see (1.2), as indeed clearly entailed by ()3.4|) with (1.5). 

The actual behavior of these similarity solutions z (t) need not be discussed any further, 
nor explicitly displayed, since the interested reader will have no difficulty to figure it out 
from (|3.1a|) (or (|3.4j) ) and (1.6). Let us only note that, in the "equal particle" case in 
which all the coupling constants a nm coincide, 

a nm = a i (3.5a) 
the expression (|3~3d|) yields 



Cn(r) 



c n exp (Pt) , 



(3.4) 



A = a(M-l), 



(3.5b) 



hence (|3.2b|) yields (provided A / -1; otherwise, see the previous paragraph) 



r = [l + a(M-l)] 



-l 



(3.5c) 



and in this special case the constants c n can be explicitly exhibited, 



c n = cexp (2inn/M) 



(3.5d) 
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consistently, see (|H.^a,j) and (|3.5b|) . with the identities 

N 

V exp(2ivrm/M)/ [exp(2i7rn/M) - exp(2zvrm/M)] = — (M - l)/2, 

m=l; mj^n 

n = l,2,...,M. (3.6) 

The treatment of this special "equal particle" case is not new [4, 2]. 

Let us finally note that, even though in general the constants c m cannot be determined 
in explicit closed form in the general case with arbitrary N and different coupling con- 
stants a nm , whenever these depend symmetrically from their two indices, see (1.12), there 
holds the sum rule 

N 

£> n = 0, (3.7) 

n=l 

which clearly obtains by multiplying (|3.3aj) by c n and summing over the index n from 1 
to N (since the double sum in the right-hand side then vanishes due to the antisymmetry 
of the summand). For N = 2 this sum rule is sufficient to determine the constants c n (up 
to a common rescaling constants), since it clearly entails c\ = c, c 2 = — c. 



4 The two-body problem 

In this section we provide a somewhat more complete treatment of the two-body problem 
than given in [1] . This is of interest in itself, but even more so for the insight it provides, not 
only for the two-body case but as well for the iV-body case (as discussed in the following 
Section 5), on the nature of the singularities of the solutions ( (t) of (1.7) as functions of 
the complex variable r, hence on the periodicity of the corresponding solutions z (t) of the 
"physical" equations of motion (1.5). 

For N = 2 the equations of motion (1.7) are consistent with the assignment (corre- 
sponding to the standard separation of the center of mass and relative motions) 

Ci (r) + C2 (r) = Ci (0) + C2 (0) + Vr, (4.1a) 
Ci(r)-C 2 (r) = C(r), (4.1b) 

namely 

Ci (t) = [Cl (0) + C2 (0) + Vr + C (r)] /2, (4.1c) 
C2 (r) = [Ci (0) + C2 (0) + Vr - C (r)] /2, (4. Id) 

where (see (1.14), (1.8) and (|4~Ta|0 

V = Ci (0) + C 2 (0) = Cl (r) + C 2 (r) = 2Z (0) (4.2) 
is a (generally complex) constant and the difference C (t) satisfies the second-order ODE 



C" = a 



V 2 -{C'f /C (4.3) 
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Here and throughout this section a = a\% = 021 is the relevant "coupling constant" , and 
primes denote of course differentiations with respect to r. 

This ODE is easily integrated once (after multiplying it by the factor 2£' /[V 2 — (C) 2 ])i 
and one gets thereby 

(C) 2 = V 2 + BC 2a , (4.4a) 

( C ') 2 = y 2 [i + (c/L)- 2a ], (4.4b) 

with B = — 4£[ (0) C2 (0) [C (0)] 2a - or, when V 7^ 0, equivalently but notationally more 

conveniently, L = (B/V 2 ) 1/{2a) = C (0) |[C' (0) /V} 2 - ^j. 1 2< ^ _ a (generally complex) 
integration constant. 

In the special case V = the first-order ODE Q4.4aJ) can be easily integrated once more, 
and one obtains thereby [1, 2] the solution of (|4.3|) in closed form: for a ^ — 1, 

C(r) = c(r-nr (4.5a) 

with 

7 = 1/(1+ a); (4.6) 
for a = —1, 

C(r) = cexp(/3r). (4.7a) 

Here c, t\, and (5 are arbitrary (complex) constants, and it is easily seen that t respectively 
(3 are related to the initial data by the relations 

r b = -(l + a)- a C(0)/C'(0) (4.5b) 
respectively 

/3 = C'(0)/C(0). (4.7b) 
Note that this solution coincides with the similarity solution (3.1) (with N = M = 2, 

7 = r). 

For the ODE (4.4) can of course be generally integrated by a further quadrature, 

but the corresponding formula in terms of the hyper geometric function F (A, B;C;Z) [6], 



CF (l/2, -1/ (2a) ; 1 - 1/ (2a) ; [C/^]" 2a ) = V (r - r b ) , (4.8) 

is not particularly enlightening, except in the special cases listed below in which the 
hypergeometric function reduces to elementary functions (actually in these cases direct 
integration of (4.4) is the neatest way to get the solution, without going through the 
hypergeometric function) . 
For a = — 1, one easily finds 



C(r) =Lsmh[(V/L) (r - r )] , 



(4.9) 
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with To a (complex) constant (related of course to the initial data: £(0) = — L sinh[(V/L)To]). 
Hence in this case ( (r) is an entire function of r, and via (1.6) this entails that all solutions 
of the Newtonian equations of motion (1.5) (with N = 2) are in this case nonsingular and 
completely periodic with period T, see (1.2) (as entailed by (4.1) with (4.7) or (|4.9j0 . 

For a = 1 (namely, for the special value of the coupling constant that corresponds, in 
the equal-particle A^-body problem, to the solvable case, see (1.17)), it is as well easy to 
get 



with the two (complex) constants r± related to the initial values by the formula 



Note the square-root branch point in the right-hand side of (|4.1Ua|) . and its consistency 
with (4.5) and ()4.6|) (of course with a = 1). It is easily seen, via (1.6), that these find- 
ings entail the following results for the solutions z n (t) of the Newtonian equations of 
motion (1.5) (with N = 2 and a = 1). If V = (namely, if the center of mass does 
not move - initially, hence as well throughout the motion, see (1.15)), the solutions are 
nonsingular and completely periodic with period T, see (1.2), if (the initial data entail, 
see (|4.5b|) and (1.8), that) r\, falls, in the complex r-plane, outside the circular disk C 
(centered at r = i/u) and of radius 1/w); they are nonsingular and completely periodic 
with period T = 2T if T& instead falls inside the circular disk C; while they are singular 
(due to the occurrence of a two-body collision) at the time t c defined mod (T) by the 
formula rj = [exp (iut c ) — 1} / (iu), if T& falls just on the boundary of C, namely on the 
circular contour C (which is just the condition necessary and sufficient to entail that t c , 
as defined above, be real). Likewise, if instead V ^ (in which case the center of mass 
moves on a circular orbit, see (1.15)), the solutions are nonsingular and completely periodic 
with period T, see (1.2), if (the initial data entail, see (|4.10b|) and (1.8), that) the two 
complex constants t± fall either both inside or both outside the circular disk C; they are 
nonsingular and completely periodic with period T = 2T if one of the two constants t± 
falls inside the circular disk C and the other outside it; while they are singular (due to 
the occurrence of a two-body collision) if r + or r_ falls just on the boundary of C, namely 
on the circular contour C. These results of course confirm those entailed by the resolvent 
formula (1.17), as described after that formula, see above. (If one considers initial data of 
type (1.10b) with fixed z n (0) = Cn (0) and z n (0) = Q' n (0) = [iv n entailing V = fi (v± + i^) 
with fixed v n , then (|4.10b|) yields t± = t± (//) = pT X T± (1), hence the motion is completely 
periodic with period T, see (1.2), both for sufficiently small and for sufficiently large values 
of the positive scaling parameter fi; it is instead completely periodic with period T = 2T 
for a finite intermediate interval of values of //). 

A third neatly solvable case obtains if a = —1/2. Then 



C(t) = V[(t-t+)(t-t. 




(4.10a) 



r± = & (0) - Ci (0)] { [C[ (0)] 1/2 ± i [(' 2 (0)] 1/2 }' /V 2 



(4.10b) 




(4.11) 



which clearly entails that Q (r) is entire (indeed, just a second-degree polynomial) in r. 
Note that the same conclusion is implied in this case by H4.5aj) with (|4.6|) . Hence we 
conclude that in this case all solutions of the Newtonian equations of motion (1.5) (can 
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be exhibited in explicit closed form and) are nonsingular and completely periodic with 
period T, see (1.2). 

The next two cases the solutions of which we report explicitly are characterized by 
a = —2 and a = —3, respectively (the alert reader can easily check these results, or derive 
them via the approach of Appendix B). 

For a = —2 the solution reads 

C (r) = AC (0) {cn ({V (r - r ) / [AC (0)]} , 2' 1 ' 2 ) . (4.12a) 

Here and below cn(u, k) is the standard Jacobian elliptic function (see for instance [6]), 
and the two constants A and tq are related to the initial data as follows: 

A = 2- 1 / 4 [(' (0) /V] ~ 1/2 , (4.12b) 
cn ({Fro/ [AC (0)]}, 2^/2) = A. (4.12c) 

These formulas entail that, in this case with N = 2 and a = a\i = = —2, all solutions 
Cn (t) of the equations of motion (1.7) are meromorphic functions of r (note that this 
is also the case if V = 0, because a = —2 entails 7 = —1; see (4.5) and (|4.6[l ). Hence, 
in this case with N = 2 and a = a\2 = 021 = —2, all nonsingular solutions z n (t) of 
the "physical" equations of motion (1.5) or (1.1) are completely periodic functions of the 
time t, with period T , see (1.2). 
For a = — 3 the solution reads 

C (r) = AC (0) {p ({V (r - r ) / [AC (0)]} ; 0, -4)}" 1 / 2 , (4.13a) 

C (r) = V [p ((r - to) ; 0, -4 [V/( (0)] 6 {-1 + [? (0) /V} 2 })] ^ . (4.13b) 

Here and below p (it; g2, 53) is the standard Weierstrass elliptic function (see for instan- 
ce [6]), and the two constants A and To are related to the initial data as follows: 

A = {-1+[C'(0)/Vf p /6 , (4.13c) 
A 2 P ({Ft /[AC(0)]};0,-4) = 1 (4.13d) 

(the equivalence of Q4.13a|) and (|4.13bl) is entailed by the rescaling properties of the Weier- 
strass function, see for instance eq. 14.13(2) of [6]). These formulas entail that, in this case 
with N = 2 and a = a\2 = a>2i = —3, the only singularities featured by the solutions C (t) 
of the equations of motion (1.7) as functions of r are branch points with exponent —1/2 
(note that this is also the case if V = 0, because a = —3 entails 7 = —1/2; see (4.5) 
and (|4.6j> ). Hence, in this case with N = 2 and a = a\2 = °2i = —3, all nonsingular 
solutions z (t) of the "physical" equations of motion (1.5) or (1.1) are completely periodic 
functions of the time t, with period T, see (1.2), or T = 2T. 

The last case we treat in detail is characterized by a = —1/4. It can then be shown 
(again, via simple calculations analogous to those carried out in Appendix B) that in this 
case the general solution of (|4.4b|) (with V 7^ 0) is given by the (implicit) formula 

[C (r) /Lf 2 - 3 [C (r) /L] = (9/16) [(V/L) (r - r b )] 2 - 3 (V/L) (r - r b ) , (4.14a) 
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entailing 



oo 



C(t) = V(t 



n) £ ai [(V/L) ( 



T — 




(4.14b) 



with 



CtQ = Ot\ = 1 



«2 = -3/16, 



(4.14c) 



It is thus seen that, in this case with a = —1/4, £ (r) has a square-root branch point at 
r = Tb if the center of mass moves (V 7^ 0), while it has a cubic-root branch point if 
the center of mass does not move (V = 0; indeed a = —1/4 entails 7 = 4/3, see (|4.5aj) 
and (14.61) ). Hence, in this case with a = —1/4, the solutions z n (t) of the physical two- 
body problem (1.5) are nonsingular and completely periodic with period T, see (1.2), if 
the branch point (whose location depends of course on the initial conditions, see l)4.5b|) 
or (|4.14a|l ) fall outside the circular disk C (centered in the complex r-plane at r = i/u 
and of radius 1/w); they are of course singular at a finite real time t = t c (when the 
two particles collide) if r& falls just on the boundary of C; while if the branch point r& 
falls inside the circular disk C they are nonsingular and completely periodic with period 
T = 2T or T = 3T depending whether the center of mass does (V 7^ 0) or does not 
(V = 0) move. 

As implied by our general treatment and exemplified by these examples, an analogous 
analysis of the emergence of periodic motions in the two-body problem (1.5) can also 
be made for an arbitrary value of the coupling constant a. The conclusions of such an 
analysis depend essentially on the nature, and location, of the branch points featured by the 
solutions of the first-order ODE (4.4) (with V 7^ 0; the V = case is completely illuminated 
by the explicit solution (4.5) with (|4.6j) . hence our treatment below refers exclusively to 
the V 7^ case). A detailed treatment of this question is given in Appendix B, and it 
leads to the following results. 

Three cases must be distinguished, depending on the value of the real part of the 
coupling constant a. 



In this case the branch point at, say, r = t^, is characterized by the exponents 7 and 1 — 7, 
see (|4.6j) . with the behavior of £ (r) for r ~ t\, given by the formula 



(for a justification of this formula, including the significance of the coefficients gtu see 
(B.20) with (B.2)). Note that this formula entails C ( r fe) = ( as wen as IC'( r fe)l = °°i 
since (|4.6j) and (j4.15aft entail < Re (7) < 1), namely this singularity is associated with a 
"collision" of the two particles (see (|4.1b|l ). which both move with infinite speed when they 
collide. This singularity is analogous to that which obtains in the V = case, see (|4.5a|) . 



Case (i): 



Re (a) > 0. 



(4.15a) 



ar) = L 1 -^{V/L){T-n)V 




(4.15b) 
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Case (it): 



1 < Re (a) < 0. 



(4.16a) 



In this case the branch point at, say, r = r&, is characterized by the exponent 



(3 = -2a, 



(4.16b) 



with the behavior of £ (r) for r ~ r& given by the formula 



C (r) = ±y (r - r b ) \l + J2 91 [± (V/i) (r - n)] 



i=i 




(4.16c) 



This singularity, however, is generally not of the same type as that which obtains in the 
V = case, see (|4.5a|) (except for those values of a such that 7 = (1 + a) -1 and = —2a 
differ by an integer). In both cases (V = 0, see (|4.5ap with (|4.6|) and (|4.16aj) : 7/0, see 
(4.16)) C ( T b) = 0, namely this singularity is again associated with a "collision" of the two 
particles (see (|4.1b|0 : however, in the V = case \(' (Tb)\ = 00, while in the V ^ case 
(' (77,) = ±V and this entails (see (4.1)) that either (r) or £2 (t) vanishes at r = so 
that 



In this case both behaviors, (|4.15b|) respectively (|4.16c|) . are possible, so both type of 
branch points, characterized by the exponents 7 = (1 + a) -1 respectively (5 = —2a, may 
be present: but the first type of branch point, characterized by the exponent 7 (which 
obtains now for all values of V) corresponds now to the phenomenon of "escape to in- 
finity" (|C( r fe)| = lC'( r fe)l = °°)> while the second type of branch point, characterized by 
the exponent (3 (which obtains only if V does not vanish, V 7^ 0) corresponds to the 
phenomenon of "two-body collision" (Q (r^) = 0, with (' (17,) = ±V entailing (|4.16d|0 . 

The behaviors of the corresponding solutions, z (i), of the "physical" equations of mo- 
tion (1.5) depends on the locations of the branch points r?, in the complex r-plane, which 
of course depend themselves on the initial data (and on the coupling constant a). (Note 
that the above analysis is local, namely it applies in the neighborhood of each branch 
point; nothing excludes that there be several, or possibly an infinity, of them). If none 
of these branch points is located inside the disk C (centered in the complex r-plane at 
r = i/uj and of radius 1/uj) nor on its boundary C, then the solutions z n (t) are nonsin- 
gular and completely periodic with period T, see (1.2) (irrespective of the value of the 
coupling constant a). If instead one branch point, say at r = t^, falls inside the disk C, 
then the motion is again nonsingular but generally not periodic, unless the branch point 
exponent is (real and) rational, in which case the motion may again be nonsingular and 
completely periodic (see below for a justification of the conditional), but with a larger 
period T = qT, where q is the (positive) integer in the denominator of the rational expo- 
nent characterizing the branch point, which depends of course on the two-body "coupling 



Ci(n)C 2 (n) = o. 



(4.16d) 



Case (Hi): 



Re (a) < -1. 



(4.17) 
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constant" , as discussed above (indeed a necessary and sufficient condition for the branch 
point exponent to be rational is that the coupling constant a be itself rational - except for 
the special case a = —1, when the motion, see (4.7) and (|4.9j) . is always nonsingular and 
completely periodic with period T, see (1.2)). And if falls just on the boundary of C, 
namely on the circular contour C, then at a finite real time t c defined mod (T) by the 
formula Tj, = [exp (iu;t c ) — 1] / (iuj) the "physical" equations of motion (1.5) become gen- 
erally singular, either because the two particles collide, or because they escape to infinity 
(as the case may be, see above). Note however that, if a is a negative integer or a negative 
half-integer, in the case in which the two particles collide there is no singularity; this is 
entailed by (|4.1(icj) . it is consistent with the occurrence of a collision thanks to (|4.16d|) . 
and it is indeed verified in the special cases a = —2, a = —3, a = —1/2, see (4.12), (4.13), 
(4.14). If more than one branch point occurs inside the disk C the analysis must be ad- 
justed accordingly; of course the outcome is critically affected not only by the presence of 
such branch points, but as well by which sheets they are located on, namely it depends on 
the overall structure of the Riemann surface associated with £ (r), the key element being 
always the path traveled on that surface by the complex point r = [exp (icut) — 1] / (iuo) 
as the real variable t ( "time" ) evolves onward from the initial moment t = 0. 

In conclusion let us emphasize that our findings identify the following reasoning as 
instrumental to understand the motion of the particles in the plane entailed, for given 
initial data, by the Newtonian equations of motion (1.1), or equivalently by the complex 
equations of motion (1.5). The idea is to fix attention on the solution of (1.7) (rather 
than (1.5)) corresponding to the same initial data (see (1.8)). This defines the solution 
Cn ( r ) (with n = 1,2, since we are now restricting attention to the two-body case; but 
we will see in the next section that essentially the same reasoning applies in the N-body 
case), to which is generally associated a multi-sheeted Riemann surface in the complex 
r-plane. The behavior of the solution of the "physical" equations of motion (1.5) as a 
function of the, of course real, time t is then obtained by traveling on that Riemann 
surface following the circular contour C defined by (1.6). Depending on the structure of 
the Riemann surface, this may entail a motion that is nonsingular and completely periodic 
(with period T, see (1.2), or with a period which is an integer multiple of T), that is 
singular (if a branch point happens to sit just on the contour C), or that is nonsingular 
but not periodic. Two mechanisms may give rise to the latter outcome (no periodicity): (i) 
the nature of the branch points, if they are characterized by an exponent that is not a real 
rational number (whether this is going or not to happen is immediately predictable, since 
it depends on whether the coupling constant a is or is not itself a real rational number, see 
(IT1T1) and (14.1 6hjH : (%%) even if the coupling constant is a real rational number, so that each 
branch point yields only a finite number of sheets, there still may be an infinite number of 
sheets due to an infinity of branch points (of course not all of them occurring necessarily 
on the same sheet, but possibly in a nested fashion), and it may then happen that by 
traveling along the contour C an endless sequence of new sheets is accessed (a necessary 
condition for this to happen is that the Riemann surface feature an infinity of branch 
points inside the contour C - of course, on different sheets). Of course both mechanisms 
could be at work simultaneously. 

Clearly the completely periodic motions correspond to an integrable behavior of the 
system; and, as shown above, whatever the value of the coupling constant a, there al- 
ways exist at least a set (having nonvanishing, indeed infinite, measure in phase space) 
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of initial data which yield such a behavior (nonsingular and completely periodic motions 
with period T, see (1.2)). As for the motions which are instead not periodic, one may 
well ask whether they should be categorized as chaotic or as not chaotic. It seems natural 
to suggest that, if the chaotic behavior is characterized by the fact that two trajectories, 
however close the initial data are that define them, give eventually rise to completely 
different motions ("sensitive dependence on initial conditions"), then clearly the nonpe- 
riodic motions will be chaotic respectively not chaotic according to whether it is infinite 
respectively finite the number of branch points which, by falling just inside or just outside 
the circular contour C, determine the (infinite) number of sheets that are accessed by a 
path following that contour C on the Riemann surface associated with the solution £ (r) 
of (1.7). Indeed, if that number of branch point is finite, two sets of initial data that 
are sufficiently close to each other (in phase space) yield two Riemann surfaces which are 
sufficiently similar to each other so that, throughout the two motions corresponding to 
the two sets of initial data, the same sequence of sheets is traveled. But if that num- 
ber of branch points is infinite, then even two solutions C ( r ) °f (1-7) that are initially 
very close will inevitably be associated with two Riemann surfaces in which one relevant 
branch points falls in one case inside, and in the other outside, the circular contour C, 
hence, by traveling on C, after that point has been passed the corresponding trajectories 
of the physical problem (1.5) (or, equivalently, (1.1)) become different, because from that 
moment a different sequence of sheets is accessed of the two Riemann surfaces associated 
with the corresponding solutions Q (r) of the evolution equations (1.7). This is then to be 
interpreted as the cause of chaos. Note however that, even in the chaotic case, there is no 
local exponential divergence of trajectories; the mechanism that causes the onset of chaos 
in this case is rather analogous to that which characterizes the nonperiodic free motion of 
a point in, say, a triangular plane billiard with angles which are irrational fractions of tt 
(then any two trajectories, however close they initially are - and for however long they 
remain close - eventually become topologically different because one of the two misses 
a reflection that the other one takes, and from that moment onwards their evolutions 
become quite different). 



5 Branch points of the solutions of (1.7) in the iV-body case 

In this Section we investigate the branch point structure in the complex r-plane of the 
solutions C (t) of the equations of motion (1.7), since - as discussed above - the nature and 
location of these branch points determine the behavior of the solutions z (t) of the equations 
of motion (1.5), namely of the Newtonian equations of motion (1.1), as functions of the real 
time variable t. Let us re-emphasize that it is indeed the structure of the Riemann surface 
associated with the solution ( (r) of the equations of motion (1.7) that determines whether 
the corresponding solution z(t) of the "physical" equations of motion (1.5) namely (1.1) 
does or does not become singular as function of the real time variable t, and if it is not 
singular throughout time whether or not it is completely periodic, and if it is periodic 
then with what period (whether T, see (1.2), or an integer multiple of T). The rule to 
evince these conclusions is quite simple, see (1.6): to obtain z(t) as function of the real 
variable t one must follow the corresponding solution ( (r) (namely, that characterized 
by the same initial data, see (1.8)) as the complex "time-like" variable r travels, on the 
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Riemann surface associated with that solution £ (r) , around and around counterclockwise 
along the circular contour C centered at i/u and of radius 

The implications of this analysis have already been discussed in the preceding Section 4 
in the context of the two-body problem. The situation in the iV-body case is somewhat 
analogous. Indeed the structure of the evolution equations (1.7) entails that the same two 
mechanisms discussed in the two-body case are generically responsible for the emergence 
of singularities, at some complex values t&, of the solutions of these equations of motion, 
(1.7), for an arbitrary number N of particles: namely, singularities arise either from the 
"collision" of two particles, characterized by the relation £i (77,) = C2 (t&) (with \d (77,) | = 
IC2 (t&)| < 00, where we assign, without loss of generality, the labels 1 and 2 to the two 
"colliding" particles) , or from the simultaneous "escape to infinity" of two or more (say, M) 
particles, characterized by the relation \Q\ (77,) = |^2 (t^)| = • • • = \Qm (t&)| = 00 (where, 
without loss of generality, we assumed the M "particles" ( n (r) that escape to infinity 
as t — > t f, to be labeled by the first M indices, n = 1, . . . , M). Note that here we use 
again inverted commas around the word "collision" and "escape to infinity" to underline 
that only in the special cases in which, via the transformation (1.6), to the value 77, there 
corresponds a real value t c (namely, 77, = [exp {iujt c ) — 1] / (iuj) with t c real; see (1.6)), the 
"collision" , or the "escape to infinity" , corresponds indeed to a real event for the physical 
problem (1.5) namely (1.1). 

The statement we just made is not meant to exclude the possibility that "collisions" 
involving simultaneously more than two particles occur and cause correspondingly a sin- 
gularity: indeed the exact similarity solutions (with M > 2) presented in Section 3 provide 
examples of solutions characterized by such phenomena. But it stands to reason - and 
it is confirmed by our analysis, see below - that M-particle "collisions" with M > 2 are 
not generic, namely they are not associated to solutions emerging from the assignment 
of generic initial data: indeed for a generic solution of the equations of motion (1.7) the 
complex equation (1 (77,) = C2 (t&) generally has at least one complex solution 77, (and more 
likely many such solutions, indeed quite possibly an infinity of them), while one should 
not expect the M — 1 complex equations ("1 (t&) = C2 (t&) = • • • = Cm ( T b)i with M > 2, 
to possess any solution at all (although there are of course special solutions £ (r) of (1.7) 
for which such multiple equations do possess solutions, see for instance Section 3). Let us 
then understand the type of singularity associated to these two types of "events" . To this 
end the discussion of the preceding Section 4 is helpful (especially to guess at the nature 
of the singularity associated with such "events"), but a new treatment in the iV-body 
context is nevertheless necessary. 

We analyze firstly the singularities associated with "two-body collisions" . 

For notational convenience let us assume, without loss of generality, that the two par- 
ticles involved in the two-body event - which happens at r = 77, - carry the labels 1 and 2, 
and let us call a the coupling constant associated with this particle pair, 

a = a±2 = a 2 i- (5.1) 

Let us firstly assume that 



Re (a) > 0, 

so that the real part of the branch-point exponent 7 (see (4.6), and below), 
7 = 1/(1 + a), 



(5.2a) 



(5.3) 
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satisfies the restriction 

0<Re(7)<l. (5.2b) 

It can then be shown (see Appendix C) that in the neighborhood of a "two-body collision" 
occurring at r = r b the solution £ (r) of (1.7) features the following behavior: 

Cs (r) = 6 + (-I)*" 1 c (r - r b y + v(r- r b ) 

oo oo 

+ £ £ 41 (— T^-M, , = 1,2, (5.4a) 

k=l l,m=0; l+m>l 

Cn (r) = &„ + u„ (r - r 6 ) 

oo oo oo 

+ r b r^^ , n = 3, . . . , N. (5.4b) 

k=l l=S kl m=0 

In these formulas the 2N (complex) constants b, c, v, r b and b n , v n (with n = 3, . . . ,N) 
are arbitrary, except for the inequalities 

b^b n , b n ^b m , n,m = 3,...,N, (5.4c) 

while the coefficients of the sums are determined in terms of these constants (see Ap- 
pendix C). Clearly these formulas, (5.4), entail Ci ( r f>) = C2 (jb) = b Cn (t&) = b n 7^ 
Cm (n) = b m for n / m, 3 < n,m < N with |Ci (n)\ = IC2 (^6)1 = H < but (see 
(j5.4aj) and (|5.2b|) ) |C( (r 6 )| = | C2 ( r fc)l = °°> while |Cn (tj,)| = \b n \ < 00 and (see (|5.4b|) and 
(|5.2b|l ) |C4 (Tfc) I = \v n \ < 00, n = 3, . . . , N. This confirms that the corresponding event is 
to be interpreted as a "collision" of the two particles 1 and 2, with both colliding particles 
moving infinitely fast at the collision time r = r b . What interests us most is the nature 
of the corresponding singularity: a branch point characterized by the exponents 7 and 
1 — 7, see ()5.3|) . And the fact that such a singularity is associated with an expression of 
the solution C ( T ) °f (1-7) that features, see (5.4), the maximal number, 2iV, of arbitrary 
constants, demonstrates the generic character of such singularities, which are therefore 
likely to be featured by the solutions ( (r) corresponding to a generic set of initial data. 

Of course a more complete treatment shall have to face the nontrivial task to prove 
that the iteration used in Appendix C to justify (5.4) does converge and thereby yields 
a solution which has indeed the analytic structure suggested by the formulas written 
above. In this paper we limit our consideration, both here and below, to this incomplete 
treatment, which is however sufficient to guess the character of the branch points. 

Likewise, if the inequality (|5.2aj) is reversed, 

Re(a)<0, (5.5a) 

it can be shown (see Appendix C) that in the neighborhood of a "two-body collision" 
the behavior of the solutions C ( r ) i s characterized, rather than by (5.4), by the following 
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expressions: 



Cl ( T ) = 6 + c (r-r 6 ) 1+/3 + Yl guir-V,)^, (5.6a) 

k,l=l;k+l>3 



oo oo 



C2 (r) = b + V2 (r - r 6 ) - c (r - r b ) 1+/3 + ^ ^ (r - t 6 )*+^ , (5.6b) 

fc=l i=25 fcl 



oo oo 

Cn (r) = 6n + (r - r 6 ) + £ £ ^ (r - r b ) fc+ ^ , n = 3, . . . , N, (5.6c) 



k=2 1=0 



with 



= -2a (5.7) 
so that (see ()5.5a|) ) 

Re (p) > 0. (5.5b) 

In these formulas the 2N (complex) constants 6, c, V2, and b n , v n (with n = 3, . . . , iV) 
are arbitrary (except, again, for the inequalities (|5.4c|0 . while the coefficients of the sums 
are determined in terms of these constants (see Appendix C). Clearly these formulas 
(see (5.6) and l)5.4cjl) entail again Ci(n) = bin) = b / ( n (n) = b n ^ C m (n) = b m 
for n ^ m, 3 < n, m < N with |£i (tj,)| = IC2 (t&)| < 00, but now (see ()5.6a|) with ()5.5b|) 
and (|5.6b|l ) (r^) = 0, (' 2 (r(,) = ^2 (so that ([ (rb) C[ i T b) = 0; see the right-hand side 
of (1.7)), while of course again |£ n (r^)! = \b n \ < 00 and (see (|5.2b|0 |C4 (^fe)l = \v n \ < 00 
for n = 3, . . . , N. Hence the corresponding event is again to be interpreted as a "collision" 
of the two particles 1 and 2, but now with particle 1 having zero velocity at the time of 
the collision, in contrast to particle 2 which moves with velocity i>2 (note the notational 
distinction thereby introduced among the two colliding particles); while for our purposes 
the interpretation of these formulas, (5.6a,b,c), is that the nature of the corresponding 
singularity is a branch point characterized by the exponent /?, see (|5.7j) . And again the fact 
that such a singularity is associated with an expression of the solution Q n (r) of (1.7) that 
features, see (5.6a,b,c), the maximal number, 2A r , of arbitrary constants, demonstrates 
the generic character of such singularities, which are therefore likely to be featured by the 
solutions C (t) corresponding to a generic set of initial data. 

This concludes our discussion of the singularities associated with "two-body collisions" . 
Those associated with "collisions" involving more than two particles can be discussed in 
an analogous manner (also taking advantage of the results of Section 3), but in view of 
their lack of genericity we forsake their treatment here, and we rather proceed to discuss 
the singularities associated with "escapes to infinity". 

This phenomenon can only occur if, for some group of the interacting particles, which 
without loss of generality is hereafter assumed to be formed by the M particles with labels 
from 1 to M (where 2 < M < N), the quantity A, 

M 

A = (2/M) £ a nm , (5.8a) 
n,m=l; n>m 
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has real part less than negative unity, 

Re(A)<-l, (5.8b) 
so that the corresponding quantity T (see (3.2b)), 

r = (l + A) _1 , (5.9a) 
has negative real part, 

Re(r)<0. (5.9b) 

The quantities a nm appearing in the right-hand side of (|5,8a|) are of course the coupling 
constants that characterize the two-body interactions acting among the particles belonging 
to this group of M particles, see (1.7) or (1.5); and there may of course be, for a given 
A^-body problem, several subgroups of particles such that the corresponding quantity T, 
defined according to the above prescription, has negative real part, see (|5.9b|) . Let us 
assume for simplicity that there is just one such group. It is then easily seen that the 
dominant term at r ~ T5 of the solution £ (r) representing the "escape to infinity" of 
the M coordinates Q n (r), n = 1, . . . , M, reads 

Cn (r) ^c n {T-r b f , n = 1, . . . , M, (5.10a) 
Cn(r) = 6„, n = M + l,...,N. (5.10b) 

In these formulas the N (complex) constants b n are arbitrary, while the coefficients c n , 
n = 1, . . . , M are determined (up to a common rescaling factor) in terms of the coupling 
constants a nm with n, m = 1, . . . , M by (3.3a) with (|5.8aj) . Clearly these formulas (see 
(5.10) and (OEll I entail [£„ (r 6 )| = \(' n (r 6 )| = 00 for n = 1, . . . , M and Cn (n) = b n for 
n = M + 1, . . . , N. Hence the corresponding "event" at r = is indeed to be interpreted 
as the "escape to infinity" of the M particles with (conveniently chosen) labels from 1 
to M. But, as already indicated in the introductory Section 1, a discussion of the analytic 
structure of this solution £ (r) at r « 17, is postponed to a separate paper. 



6 Analysis of various motions 

In this section we survey the implications of the findings discussed above (especially in 
Section 5) for the motions of the physical models (1.1) characterized by various choices 
of the number of particles and of the coupling constants a nm . We always refer, for 
notational convenience, to the complex avatar (1.5) of these equations of motion, and we 
display the trajectories of solutions produced via a computer code manufactured by one 
of us (MS) for the numerical integration of these equations of motion, (1.5). Additional 
numerical examples, as well as the discussion and display of the numerical code, will be 
reported elsewhere [8], enabling thereby any interested party to enjoy an ampler personal 
exploration of analogous numerical results (easily generated on a personal computer). 

The basic mode of presentation of our results begins by choosing a specific model, 
characterized by the number A^ of particles (for simplicity we restrict consideration below 
to three- and four-body cases, N = 3 and N = 4, although the simulation program can 
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actually handle an arbitrary number N of moving points) and by a set of coupling con- 
stants a nm , see (1.5), which we always assume to satisfy the symmetry property (1.12). We 
moreover always set, for notational convenience (and of course without loss of generality, 
see the Remark II .1(1 

to = 2ir, (6.1a) 

entailing an assignment of the time scale such that the basic period be unity, 

T=l, (6.1b) 

see (1.2). We then generally consider, for each model, a sequence of solutions of the 
equations of motions (1.5) corresponding to a sequence of initial data characterized by 
the same initial configuration of the particle positions in the plane and by sets of initial 
velocities modified, from one solution to the next one considered, via multiplication by an 
increasing sequence of positive factors [i (by analogy with the notation in (1.10b)): 

z n (0) = 4°\ i n (0)=/i4°). (6.2) 

As implied by the treatment of the preceding sections (see in particular Section 1 and 
the Remarks 11.11 and I1.3JI this entails that all the sequentially analyzed solutions z (t) of 
the "physical" equations of motion (in their avatar (1.5), and characterized by a specific 
assignment of iV and a nm ) correspond, via the fundamental relation (1.6) (with an appro- 
priately chosen value of the constant u in this formula (1.6), namely u = 2-/r//i), to the 
same solution ( (r) of (1.7), characterized by the same values of N and a nm , and by the 
same initial data, 

Cn(0) = 4 0) , C;(0) = 4 0) ; (6.3) 

hence the solutions z[t) of the equations of motion (1.5) discussed below are obtained 
via (1.6) from the same solution ( (r) of (1.7) by traveling on the Riemann surface asso- 
ciated with this solution of (1.7) on circular contours C = C (/i), in the complex r-plane, 
the diameters of which always lie on the upper imaginary axis of the complex r-plane, 
with their lower ends fixed at the origin, r = 0, while the size of these diameters is n/n, 
namely it gets increased via multiplication by the same rescaling factor, fi, that yielded 
by multiplication the set of initial velocities. The sequence of solutions we consider starts 
generally with a completely periodic one having the basic period T = 1; note that, as 
implied by Proposition ^21 an d as of course confirmed by the numerical simulations, such 
solutions z_ (t) always exist for sufficiently small initial velocities, namely sufficiently small 
values of /i (they correspond to solutions £ (r) of (1.7) which are singularity free, on the 
main sheet of their Riemann surface, inside the circle C (//) - which, for very small fj,, 
becomes very small, and gets very close to the origin, r = 0, where the solution £ (r) is of 
course holomorphic) . The sequence of solutions z{t) identified according to the procedure 
indicated above corresponds therefore generally to the sequential coming into play of dif- 
ferent branch points (see Section 5) of the same solution £ (r) of (1.7), as they get enclosed 
inside the circular contours C = C (fi) when the size of these contours gets enlarged pro- 
portionally to the parameter /j, (which is sequentially increased from one simulation to the 
next one). The inclusion of such branch points may entail the transition - for models with 
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real rational coupling constants, to which our consideration is mainly restricted below, as 
they display the most interesting phenomenology (aside from the fact that, in numerical 
computations, all numbers are rational!) - to solutions which are again completely peri- 
odic but with larger periods, or which possibly are no more periodic (although of course 
it is hardly possible to distinguish numerically a nonperiodic solution from one which is 
periodic but with a very large period, see below), according to a phenomenology which 
can be understood on the basis of the treatment of the preceding Sections 4 and 5 and 
the discussion at the end of Section 1 - an analysis which turns out to be indeed vali- 
dated by the numerical results reported below (and see also the analogous discussion in 
Ref. [9]). We also present below some results with complex coupling constants, and tersely 
discuss the associated phenomenology. On the other hand, as already mentioned above, 
our consideration is restricted to coupling constants the real part of which is positive, 
Re (a nm ) > 0. 

Enough of preliminary remarks: let us proceed to display and comment a few represen- 
tative examples. In the graphs displayed below, particle 1 will be shown in Red, particle 2 
in Green, particle 3 in Blue and particle 4 (if present) in Yellow. Whenever we felt such an 
additional indication might be usefully displayed we indicated with a black diamond the 
initial position of each particle (at t = 0), with a black dot the position at a subsequent 
time t = ti, and with smaller black dots the position at every subsequent integer multiple 
of t\ (namely at t = = kt\, k = 2,3, . . .); in this manner the direction of the motion 
along the trajectories can be inferred (from the relative positions along the trajectories of 
the diamond and the larger dot), as well as some indication of the positions of the particles 
over time, as they move (by counting the dots along the trajectory). Of course a much 
more satisfactory visualization of the behavior of the many-body system is provided by 
simulations in which the particle motions are displayed as they unfold over time (as in a 
movie); it is planned to make available soon, via the web, the numerical code suitable to 
perform such simulations on personal computers [8]. Let us emphasize that such simula- 
tions are particularly stunning to watch in the case of high-period trajectories, which are 
very complicated (see below), so that the fact that the particles return eventually exactly 
on their tracks appears quite miraculous and is indeed a remarkable proof of the reliability 
of the numerics. 

The first example we report is characterized by the following parameters: 
N = 3; au = a 2 i = a RG = 1, ai 3 = a 3 i = a RB = 1/2, 

fl23 = ^32 = a GB = 3/2, (6.4a) 

and by the following values of the parameters x$ , Vn^ , Xn^ , Vn^ characterizing the initial 
data via (|6.2|) : 



x 



(0) 



= 




= -1 



V { ? = -0.5; 



i 



x. 



(0) 



= 1 



= 0.5 




2 



X 



(0) 



= 2 



= 0.5 



(6.4b) 



3 



Our numerical findings are synthesized in the Tables 6.1a,b, and in the sequence of 
Figures 6.1 (the contents of which are identified in Table 6.1a). Here and below we employ 
the acronym HSL (Hie Sunt Leones) to denote a (presumably) chaotic behavior. We trust 
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the significance of this material - as presented here and in the rest of this section - to be 
self-evident, but we nevertheless add a few comments. 



Table 6.1a 



A* 


0.8 


1 


1.301 


1.302 


5.28 


5.29 


5.93 


5.94 


6.82 


6.83 


50 


Period 


1 


1 


1 


2 


2 


6 


6 


7 


7 


HSL 


HSL 


Fig. 6.1 


a 




b 


c 


d 


e 


f 


g 


h 


i 





Table 6.1b 



Period Change 


1 -> 2 


2^6 


6^7 


7 -> HSL 


Colliding Particles 


R-G 


G-B 


R-G 


R B 


Collision Time 


0.80 


0.77 


2.87 


3.72 



Clearly for /x < 1.301 the solution £ (r) of (1.7) (with (|6,3() and ()6.4bjl ) is holomorphic 
inside the circle C (//), hence the (corresponding) solution z{t) of (1.5) (with (6.1), Q6.2|) 
and (6.3b)) is completely periodic with period T = 1 (see the first 3 columns of Table 6.1a, 
and Figs. 6.1a,b). For 1.302 < /x < 5.28 this (same) solution £(r) of (1.7) (with (JOJ) and 
(|6.4b|l ) has a single singularity inside the circle C (/i), a branch point with exponent 



7-RG = V (org + 1) = 1/2; 



(6.5a) 



hence by traveling on the Riemann surface associated with £ (r) on the circle C (/x) two 
sheets are now visited before returning to the point of departure, and therefore the (cor- 
responding) solution z(t) of (1.5) (with (6.1), (|6.2|) and (6.3b)) is completely periodic with 
period T = 2 (see the fourth and fifth columns of Table 6.1a, and Figs. 6.1c,d). The 
transition from the regime with period T = 1 to that with period T = 2 occurs due to the 
entrance of the relevant branch point from outside to inside the circle C (/x); this corre- 
sponds to the occurrence of a collision among the Red and Green particles (see Figs. 6.1b 
and 6.1c, and the first column of Table 6.1b), consistently - according to the analysis of 
Section 5 - with the presence of the coupling constant a^c = 1 in (|6.5a|) . 

For 5.29 < /x < 5.93 this (same) solution C(r) of (1.7) (with (JOJ and (|Ohjn has two 
singularities inside the circle C (/x), namely the branch point with the exponent (|6.5aj) and 
a second one (located either on the main sheet, or possibly on the second sheet accessed 
via the cut associated with the first branch point, with exponent (|6.5a|) L with exponent 



1GB = V (a G B + 1) = 2/5, 



(6.5b) 



which opens the way to four additional sheets; hence by traveling on the Riemann surface 
associated with £ (r) on the circle C (/x) altogether six sheets are now visited before re- 
turning to the point of departure, and therefore the (corresponding) solution z{i) of (1.5) 
(with (6.1), ()6.2|) and (6.3b)) is completely periodic with period T = 6 (see the sixth and 
seventh columns of Table 6.1a, and Figs. 6.1e,f). The transition from the regime with 
period T = 2 to that with period T = 6 occurs due to the entrance (in one of the two 
sheets of the Riemann surface) of the (second) branch point, characterized by the expo- 
nent l|6.5b|l . from outside to inside the circle C (/x); this corresponds to the occurrence of 
a collision among the Green and Blue particles (see Figs. 6. Id and 6.1e, and the second 
column of Table 6.1b), consistently - again according to the analysis of Section 5 - with 
the presence of the coupling constant ogb = 3/2 in (|6,5b|) . 
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Another transition occurs for a value of [i between 5.93 and 5.94 (see the seventh and 
eighth columns of Table 6.1, and Figs. 6.1f,g), associated to another collision among the 
Red and Green particles (see the third column of Table 6.2) and bringing inside the circle 
C (fi) another branch point of square-root type (see Q6.5ajl ). that opens the way to one 
additional sheet of the Riemann surface, and thereby causes the solutions z(t) of (1.5) 
(with (6.1), (|6.2|) and (6.3b)) to be completely periodic with period T = 7 (see the eighth 
and ninth columns of Table 6.1a, and Figs. 6.1g,h) for 5.94 < fi < 6.82. 

Finally another collision occurs, for /i between 6.82 and 6.83 (see the ninth and tenth 
column of Table 6.1a, and Figs. 6.1h,i), among the Blue and Red particles (see the last 
column of Table 6.1b), which brings (on one of the seven relevant sheets of the Riemann 
surface associated with the solution ( (r) of (1.7) with (|6.3|) and (|6.4b[0 a branch point 
with exponent 



lRB = 1/ (1 + a RB ) = 2/3 (6.5c) 

inside the circle C (fj,). This opens the way to two additional sheets of the Riemann surface 
(associated with the solution ( (r) of (1.7) with (|6.3|) and (|6,4bj0 : but clearly on (at least) 
one of these two sheets there is (at least) one additional branch point (characterized by 
one of the three exponents (6.5a,b,c)) which is already located inside the circle C (/i), 
and this opens the way to additional sheets on (some of) which there also are additional 
branch points inside the circle C (fi), and so on and on, causing presumably a transition to 
chaos (see the last two columns of Table 6.1a, and Fig. 6. li) according to the mechanism 
outlined at the end of the introductory Section 1 and discussed in Section 4 (and also 
discussed, in an analogous context, in Ref. [9]). 

The second example we report is much richer. It is characterized by the following 
parameters: 



N = 3; ai2 = a 2 i = a RG = 1, a i3 = a 3 i = a RB = 2, 

a23 = 032 = a GB = 3, (6.6a) 

and by the following values of the parameters , Un^ , 4°^ , Vn^ characterizing the initial 
data via (16.21): 



*i 0) =0, y™ = t i< > = -l, #) = 1; 

*( 0) = o, vi 0) = i, 4 0) = i, ^ 0) = 0; 

4°) = -1, yf } = 0, 4 0) = -°- 5 > ^ = -°- 5 - ( 6 - 6b ) 

Our numerical findings are now synthesized in the Tables 6.2a,b, and in the sequence of 
Figs. 6.2 (the contents of which are identified in Table 6.2a). We urge the alert reader to 
repeat the analysis we detailed in the preceding example: the fact that the branch point 
exponents take in this case the values ^ R c = 1/2, ^ BR = 1/3 and job = 1/4 is of course 
relevant, and clearly some transitions from a periodic regime to another periodic regime 
with higher period are the effect of a single identifiable branch point, while others are due 
to the synergistic effect of several branch points coming into play simultaneously. The 
presence of trajectories with stunningly large periods is of course noteworthy. 
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Table 6.2a 





0.5 


0.780 


0.781 


1 


1.213 


1.214 


1.219 


1.220 


1.241 


Period 


1 


1 


6 


6 


6 


49 


49 


50 


50 


Fig. 6.2 


a 


b 


c 




d 


e 











1.242 


1.25 


1.293 


1.294 


1.295 


1.400 


1.401 


1.41 


1.442 


Period 


51 


51 


51 


57 


58 


58 


59 


59 


59 


Fig. 6.2 




f 




g 








h 







1.443 


1.721 


1.722 


1.8 


1.944 


1.945 


2.053 


2.054 


2.1 


Period 


65 


65 


66 


66 


66 


68 


68 


71 


71 


Fig. 6.2 








i 










j 





2.108 


2.109 


2.164 


2.165 


2.17 


2.171 


2.172 


2.2 


2.5 


Period 


71 


72 


72 


74 


74 


74 


HSL 


HSL 


HSL 


Fig. 6.2 










k 






m 





Table 6.2b 



Period Change 


1 -► 6 


6 ^ 49 


49 -> 50 


50 -> 51 


51 -> 57 


Colliding Particles 


R B 


R-G 


R-G 


R-G 


R-G 


Collision Time 


0.48 


3.54 


24.72 


20.59 


43.75 



Period Change 


57 -> 58 


58 -> 59 


59 -> 65 


65 — > 66 


66 -> 68 


Colliding Particles 


R-G 


R-G 


R-G 


R-G 


R-B 


Collision Time 


31.49 


27.37 


8.35 


44.25 


54.41 



Period Change 


68 -> 71 


71 -> 72 


72 -> 74 


74 -> HSL 


Colliding Particles 


G-B 


R-G 


R-B 


R-G 


Collision Time 


56.26 


5.63 


13.75 


24.905 



The third example we report is the last one with coupling constants which are both real 
and rational - namely the last one belonging to the class that can give rise to completely 
periodic motions with periods larger than unity; and it features four particles. It is 
characterized by the following parameters: 



N = 4; ai2 = C21 = org = 2, ai3 = a 3 i = orb = 5/2, 
ai4 = a.41 = ory = 1, a23 = = acB = 1/2, a24 = a42 
034 = a43 = agy = 3/2, 



acY = 3, 



(6.7a) 



and by the following values of the parameters xffl , Vt? , Xn^ , Vn^ characterizing the initial 
data via (|6.2|) : 



i 

.(0) 

; 2 

.(0) 
; 3 

.(0) 



0, 
0, 



-1.5, 



(o) 

(0) 

y 2 



l, 

0: 



(0) 
2/3 



0, 



(0) 
2/4 



0, 



x i 

~(0) 



jo) 



.(0) 

I/i 



-i; 



,(0) 



-0.5, 
= 0, 



.(0) 
2/ 2 

.(0) 

y 3 



-0.5; 



.(0) 
1/4 



(6.7b) 



The corresponding numerical results are synthesized in the Tables 6.3a,b, and in the 
sequence of Figs. 6.3 (the contents of which are identified in Table 6.3a). We again urge 
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the alert reader to repeat the analysis we detailed in the first example, taking of course 
into account that, corresponding to the six different coupling constants (see (|6.7a[) ). there 
are in this case six different branch point exponents: 7_rg = 1/3, -jbr = 2/7, jby = 1/2, 
1GB = 2/3, 7cy = 1/4 and ^by = 2/5. And the presence of trajectories with remarkably 
large periods should again be noted. 



Table 6.3a 



A* 


0.4 


0.626 


0.627 


0.692 


0.693 


1 


1.037 


1.038 


Period 


1 


1 


4 


4 


6 


6 


6 


22 


Fig. 6.3 


a 


b 


c 


d 


e 




f 


g 





1.156 


1.157 


1.207 


1.208 


1.3 


1.318 


1.319 


1.32 


1.4 


Period 


22 


23 


23 


26 


26 


26 


HSL 


HSL 


HSL 


Fig. 6.3 










h 






i 





Table 6.3b 



Period Change 


1 -> 4 


4^6 


6 -> 22 


22 -> 23 


23 -> 26 


26 -> HSL 


Colliding Particles 


G-Y 


R-G 


R-B 


R-Y 


G-Y 


B-Y 


Collision Time 


0.44 


3.33 


3.33 


9.26 


12.36 


14.53 



Finally let us display two examples with complex coupling constants, which can both be 
considered deformations of the first example considered above, obtained by adding to the 
three coupling constants characterizing that model, see (|6.4a|) . (not too large) imaginary 
parts, either all of them negative (first example below), or all of them positive (second 
example below). We also take, in both cases, the same assignment (|6.4b|) of the parameters 
that determine the initial data via (|6.2|l . 

The first of these last two examples is characterized by the parameters 

N = 3; ai2 = a 2 i = a RG = 1 - 0.3i, 

«13 = «3i = clrb = 0.5 - OAi, a 2 3 = a 32 = a G B = 1-5 - 0.2i (6.8) 

and the corresponding trajectories are displayed in the sequence of Figs. 6.4a-f; in each 
figure the value is now indicated of the scaling factor \i that determines via (|6.2|) with 
(|6.4b|) the initial data for the displayed motions, as well as the period - which is of course 
always T = 1 - in the case of periodic motions (see Figs. 6.4a,b) or instead the final time t 
of the displayed trajectories in the case of nonperiodic motions (see Figs. 6.4c-f). 

The qualitative understanding of these results is easy, and it displays an interesting 
phenomenology - limit cycles, with (asymptotic) period T = 1 - that is presumably 
characteristic of all models of type (1.5) with complex coupling constants the imaginary 
parts of which are negative, Im (a nm ) < 0. Indeed the (first) two Figs. 6.4a,b display the 
completely periodic motions, with period T = 1, which as we know are always featured 
by all models of type (1.5), as guaranteed by Proposition 11.21 and let us recall that our 
understanding of the existence of these completely periodic solutions z{t) of (1.5) with 
period T = 1 is in terms of the corresponding solution ( (t) of (1.7) being singularity 
free inside the corresponding circles C (fi) in the complex r-plane. We also know that, 
by increasing fi - hence, by increasing the size of the circle C(/x), of diameter [i/tv - 
branch points (of the corresponding solution ( (r) of (1.7)) are eventually enclosed inside 
it. These branch points are characterized by exponents 7„ m = 1/ (l + o nm ), hence by 
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exponents that, for this example, have a positive imaginary part, Im (^ nm ) > 0, since the 
imaginary part of the coupling constant a nm is negative, Im (a nm ) < 0. In this case the 
main effect of such a branch point, when one travels over the circle C (//) namely follows 
the evolution of £ (r) with r given by the formula r = [i [exp (2irit) — 1] / (2ni) (see (1.6) 
and (6.1)), is to yield the shrinking factor exp [— 27rlm {j nm ) t], as entailed by the formula 

(r - n) lnm = {fi [exp (2mt) - 1] / (2tt») - T 6 p m = exp [-27rlm ( 7nm ) t] 

x exp [2tti Re (7 nm ) i] {m [1 ~~ ex P ( — 2?ri£)] / (2ni) — tj, exp (— 27rit)} . (6-9) 

The effect of such shrinking, as long as it is relevant, is to reduce - qualitatively - by 
the same amount the magnitude of the (nonlinear) right-hand side of the equations of 
motion (1.5) or (1.7); but we know (see Section 1) that when this happens eventually 
one reaches the regime characterized by the completely periodic solutions of period T = 1 
of (1.5), the existence of which is guaranteed by Proposition 11.21 This is what one sees 
in the trajectories displayed in the Figs. (6.4c-f): in fact, given the relative scale of the 
motions of the various particles, what one apparently sees is that one of the three particles 
eventually stops (asymptotically; or it might actually tend to a periodic trajectory of very 
small size, moving along it very slowly since the period of that asymptotic motion is 
always T = 1; in either case, because of its slow motion, it essentially decouples hence 
moves on a circular orbit, neither feeling the presence of the other particles nor affecting 
their motion, due to the velocity-dependence of the two-body forces, see the right-hand 
side of (1.5)), while the other two particles approach asymptotically (in fact, given the 
graphical limitations, seem to eventually follow) periodic trajectories, always of course 
with period T = 1. Specifically, the first transition from the completely periodic regime 
to that affected by the presence of a branch point occurs at a value of /j between 1.682 
and 1.683, due to a collision among the Red and Green particles occurring approximately 
at t = 0.85 (see Figs. 6.4b,c; note that the latter of these two figures displays a motion 
in which the Green and Blue particles follow asymptotically periodic trajectories, while 
the Red particle seems to stop asymptotically). Several other transitions occur for larger 
values of /U. A representative instance, due to a collision among the Green and Blue 
particles occurring approximately at t = 0.85, is displayed by Figs. 6.4d,e: note that 
Fig. 6.4d displays a motion in which asymptotically all three particles move on periodic 
trajectories, but the Red one moves on a circle of very small size; while Fig. 6.4e displays 
as well a motion with all three particles asymptotically on periodic trajectories, but now 
with the Blue one approaching asymptotically a circular trajectory of very small size, in 
fact so small that it appears merely as a point in the figure. Finally Fig. 6.4f displays 
a motion in which the Blue particle seems to stop asymptotically, while the other two 
approach periodic trajectories. 

The interpretation of this phenomenology as outlined above is confirmed by the numer- 
ical results of the last example we report, in which the imaginary parts of the coupling 
constants a nm are instead positive, Im (a nm ) > 0. This example is characterized by the 
parameters 

N = 3; a\2 = a 2 i = a RG = 1 + OAi, 

013 = a 3 i = a RB = 0.5 + 0.3z, a 2 3 = a 32 = a GB = 1.5 + 0.2i, (6.10) 

and the corresponding trajectories are displayed in the sequence of Figs. 6.5a-f; as above, 
in each figure the corresponding value is indicated of the scaling factor /j that determines 
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the initial data via (|fi.2|) with (|6.4h|) . as well as the period - which is of course always 
T = 1 — in the case of periodic motions (see Figs. 6.5a,b) or instead the final time t 
of the displayed trajectories in the case of nonperiodic motions (see Figs. 6.5c-f). As 
in the preceding example, the completely periodic motions of period T = 1 displayed in 
the first two figures correspond to the solutions predicted by Proposition 11.21 namely 
they correspond to the solution ( (r) of (1.7) being singularity-free inside C (fi); while the 
behavior (spiraling outwardly) displayed by the other four figures is consistent with the 
coming into play of an expanding factor, the origin of which can be explained exactly as 
above (see (|6.9)) ). the difference among expanding and shrinking behavior being of course 
consistent with the change of sign of the imaginary parts of the coupling constants a nm . 
Specifically one sees from Figs. 6.5b,c that the first transition from the completely periodic 
regime is due to a collision among the Red and Green particles (occurring approximately at 
t = 0.74), and another transition is apparent from Figs. 6.5d,e, due a collision among the 
Green and Blue particles (occurring approximately at t = 0, 69). Several other transitions 
occur between these, as well as for values of \x between 1.010 (see Fig. 6.5c) and 3.474 
(see Fig. 6.5d), between 3.475 (see Fig. 6.5e) and 20 (see Fig. 6.5f), and beyond 20. 
Note that in all these simulations - with \i > 1.010, namely when the motions are no 
more periodic - all three particles eventually spiral out. This, however, does not seem 
to represent a universal behavior; other (three-body) cases have been observed, with all 
coupling constants having positive imaginary parts, in which only two particles spiral out 
asymptotically, while the third one seems to stop. Clearly a more complete analysis of this 
phenomenology will have to be made; indeed while the discussion of the last two examples 
seems to provide a satisfactory qualitative understanding of the numerical results reported, 
it also confirms the richness of behaviors displayed by the many-body problem in the plane 
(1.1), or equivalently (1.5), leaving an ample scope for additional study (for instance, of 
cases with coupling constant some of which have positive, and others negative, imaginary 
parts) - and of course there is moreover the case with coupling constants having a negative 
real part, the investigation of which shall be reported separately. 

Let us finally emphasize that the findings reported above are merely a tiny subset - 
selected on the basis of their representative character - of the data we collected in our 
numerical explorations. Indeed we expect that every interested reader will eventually like 
to engage in personal experimentations, as soon as the computer program [8] which yielded 
these results shall become available (it is planned to put it on the web soon). In particular 
we recommend once again as a thrilling experience the observation as they evolve over time 
of periodic motions with large periods, such as some of those which have been displayed 
above (but only in the guise of trajectories: see for instance Figs. 6.1e-h and especially 
Figs. 6.2c-k, 6.3c-h). 



Periodic Solutions of a Many-Rotator Problem in the Plane. II 



187 




188 F Calogero, J-P Frangoise and M Sommacal 




Periodic Solutions of a Many-Rotator Problem in the Plane. II 



189 





Fig. 6.4e 



Fig. 6.4f 
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7 Outlook 



The rich phenomenology featured by the solutions of the many-body problem in the plane 
investigated in this paper led us to refer to it, at least colloquially among us, as a "gold- 
fish" , thereby extending to this model the terminology that was initially suggested for its 
integrable variant [12]. Much additional fun and insight can certainly be gleaned from 
a more complete exploration of its solutions than we were able to report here; let us 
re-emphasize that the most illuminating way to pursue this study shall be provided by 
the computer code that simulates - by solving numerically the relevant equations of mo- 
tions, see (1.1) - the actual behavior of the many-body problem as it unfolds over time; 
this computer code is indeed a remarkably efficient research tool which we hope to make 
available soon to the scientific community at large [8]. As for the more theoretical aspect 
of the investigation of this many-body problem, let us recall that in this paper a rather 
satisfactory understanding of the behavior of the system has been provided for the case in 
which the coupling constants a nm are all nonnegative, a nm > 0, see (1.1) (or, equivalently, 
Re(a ram ) > 0, see (1.5)), while the case when this restriction does not hold still needs to 
be fully investigated. 
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Appendix A 

In this Appendix we provide a derivation of (1.9). Our starting point is the standard 
theorem that guarantees the analyticity of the solutions of (systems of) analytic ODEs 
(see for instance Section 13.21 of [7]). Our treatment is closely analogous to that of [10]. 
The standard formulation of this theorem refers to a system of first-order ODEs, say 



Wj = fj (wi,W 2 , ■ ■ .,W 2 n) , 



J = 1,2, ... , 



2 AT, 



(A.l) 



which coincides with our system, see (1.7), via the assignments 



W n = W n (r) = (n (t) - Cn (0) , 

w N+n = w N+n (r) = c [^ (r) - C (0)] , 



(A.2a) 
(A.2b) 



as well as 



fn = C 1 W N+n + C' n (0) , 



(A.3a) 



N 



[w n -W m + Cn(0) - Cm(0)] 




(A.3b) 



Here of course (as in the rest of our paper) the index n takes all integer values in the range 
1 < n < N, and N is the total number of particles. Note that with these assignments 
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all the 2N quantities Wj (r) vanish at the origin, Wj (0) = 0, as appropriate to reproduce 
the notation of [7]. In these formulas c is an arbitrary positive constant (which we have 
introduced for convenience, see below). 

The standard result [7] provides then the following lower bound for the radius p of 
the circular disk D, centered at the origin r = in the complex r-plane, within which 
the solutions w n (r) of (A.l), hence (see (A. 2a)) as well the solutions ( n (r) of (1.7), are 
holomorphic: 

p > b/ [(27V + 1) M] (A.4) 

(this formula coincides with the last equation of Section 13.21 of [7], with the assignments 
m = 2N and a = oo, the first of which is justified by the fact that (A.l) features 2N 
equations, the second of which is justified by the autonomous character of our equations 
of motion, see (A.l) with (A. 3)). The two positive quantities b and M in this equation, 
(A.4), are defined as follows. The quantity b is required to guarantee that the right-hand 
sides of (A.l) are holomorphic provided there hold the 2N inequalities 

\wj\ < b; (A.5) 

this entails, in our case, the single restriction 

b < C/2 (A.6) 

with £ defined by (1.9b), because the only source of singularities is in our case the vanishing 
of a denominator in the right-hand side of (A. 3b) (which is clearly excluded by (A.5) with 
(A.6) and (1.9b)). The second quantity in (A.4), M = M (b), is the upper bound of the 
right-hand sides of (A.l) when the quantities Wj satisfy the restriction (A.5) with (A.6); 
but of course the inequality (A.4) holds a fortiori if we overestimate M (as we shall now 
do). 

Indeed, the right-hand sides of (A. 3) clearly entail, via (A. 2) and (1.9b,c), 



M < max 



-H + 2cNa [c'H + C') 2 / (C - 2ft) 



(A.7) 



where the maximum should be taken among the two terms, while it remains our privilege 
to choose the values of the positive quantity c, as well as the value of the positive quantity b 
(respecting of course the inequality (A.6)). The positive quantity a in this formula is of 
course the maximal value of the moduli of the "coupling constants" a nm , 

a = max \a nm \ . (A. 8) 

n,m=l,...,N; n^ra 

We now make the convenient choice 

c= \(-2b(l + Na)] I \2Na('] , (A.9) 

which equalizes the two terms inside the maximum function in (A.7). Note that the re- 
quirement that c be positive entails that (A.6) must now be replaced by the more stringent 
inequality 

b < (/ [2 (1 + Na)} . (A.10) 
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Hence from (A. 7) we now get 

M <C'(t-2b)/[C-2b(l + Na) , 
and via (A. 4) there thus obtains the formula 

p > RC/C' 

with 

R = x [1 - 2 (1 + Na) x] I [(2N + 1) (1 - 2x)\ 



(A.ll) 



(A.12) 



(A.13) 



where the choice of the positive quantity x = b/Q remains our privilege, provided it falls 
in the interval 



< x < [2(1 + Na)} 



-i 



(A.14) 



which corresponds to (A. 10). 

Since (A.12) coincides with (1.9a), this latter equation is now proven. It would be easy 
to obtain the optimal choice of the constant R by maximizing the right-hand side of (A.13) 
in the interval (A.14). More simply we just take x in the middle of that interval (A.14), 
namely we set 



x= [4(1 + Na)}' 1 , 
and we thereby obtain for R the following simple expression: 
R= [4(l + 2iV)(l + 2iVa)r 1 . 

Appendix B 



(A.15) 



(A.16) 



In this appendix we investigate the nature of the branch points of the solutions of the first 
order ODE (4.4b), 



(C') 2 = V 2 \l + (C/L) 



-2a 



(B.la) 

with V and L arbitrary (nonvanishing) constants, obviously related to the initial data: 

V = c( (0) + C2 (0) (B.lb) 
(see (4.2)), and 

l/(2o) 



L = C(o){[C'(o)Af -1} 

(see (B.la)). 

Let us begin by setting 



(B.lc) 



C (r) = Lf (w) , w= (V/L) r, 



(B.2) 
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so that the equation satisfied by / (w) read 

[f'(w)] 2 = l + [f(w)]- 2a = l + [f(w)f. 

Here and below 
d = -2a, 



(B.3a) 



(B.3b) 



and the prime appended to the function / signifies of course differentiation with respect 
to its argument w. 

Two possible mechanisms are at play to yield a singularity of f (w), say at w = Wb- 
(i) the vanishing of / (w), 



f H) = 0; 
(ii) the divergence of / (w), 
f (w b ) = oo. 



(B.4a) 



(B.4b) 



Let us discuss them separately. But in both cases our approach will be based on replacing 
(B.3a) by a recursive set of equations, namely 



O'+i) f 



f H 



1 + 



Cj) 



fW 



-2a 



= 1 + 



U) 



/H 



j = 0,1,2,... 



(B.5) 



and by assigning appropriately the initial element, ^f(w), of this iteration, so that it 
account properly for the leading behavior of the solution / (w) in the neighborhood of the 
singular point w = Wb - as demonstrated by the preservation of this singular behavior 
through the iteration (see below). 

Let us begin with mechanism (i), by therefore setting 



{j) f(w b )=0, j = 0,1,2,.... 
Firstly we consider the case characterized by the inequality 
Re (a) > 0, 

and we introduce correspondingly the exponent (see below) 

7 = 1/(1 + a), 
which clearly then satisfies the restriction 

< Re (7) < 1. 
Correspondingly we set (consistently with (B.6)) 

(0) /(u;)=7 -7 ( W - Wb )\ 

and via (B.5) we then get 



(i) f 



(0) 



/H 



-2a 



1 + 



(0) 



/H 



-2a 



(B.6) 
(B.7a) 
(B.8a) 
(B.7b) 

(B.9) 

(B.lOa) 
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namely 



(0) 



/H 



1 + 



(0) 



/H 



- 2 f 2 , 



(B.lOb) 



namely (via (B.9)) 



1/2 



Wf ( w ) = y-7 ( w _ Wb y-i |i + 7 2(7-i) ( w _ w^ 1 -^ 

To get this expression we used the relation 
a 7 = 1 — 7 



which is clearly entailed by (B.8a). 

In the neighborhood of w = Wb we can rewrite (B.lOc) as follows: 



(B.lOc) 



(B.8b) 



1 + E ( 1 / 2 ) ^ (7_1) (™ - ™*f (1_7) 



i=i 



(B.lOd) 



(1) /'H = 7 1-7 (^-^) 7_1 

L/ 
I 

Here and below we use of course the usual notation for the binomial coefficient, so that 

( l [ 2 ) =r(3/2)[r(3/2-/)r(i + or 1 

= (-2)-^ (I - 1/2) (21 - 1)!!//!. (B.ll) 
We then integrate (B.lOd) (using (B.6) and (B.7b)) and we thereby get (using (B.8a)) 

^ f(w) = 7 ~ 7 (u> — Wb)" 1 

xjl + ^l + Ma)-^ l [ 2 ) 7 2,(7_1) («'-«'6) 2,(1 - 7) |- (B.12) 

A comparison of this expression, (B.12), with (B.9) demonstrates the preservation of the 
leading term characterizing the behavior in the neighborhood of the singularity. 
This result, as well as the following development, suggest setting 



(J) 



C oo I \ 

f{w) = i-\w - WbV 1 + E E - w b P(w - Wb) 21 ^ 

I 1=1 k=0 J 



j = 0,l,2..., 



of course with 



(0 W = 0, «y« = * M (l + 2Za)- 1 ( 1 / 2 ) 7 2,( ' y - 1) - 



(B.13) 



(B.14) 
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We now show that this ansatz, (B.13), is preserved by the iteration (B.5), which we 
now rewrite in the following form (analogous to (B.lOb)): 



0) 



1 + 



/H 



-2a 



1/2 



(B.15) 



w - w b ) k ~<{w - w b ) 2l{1 -^ } , 



Indeed the insertion of (B.15) in (B.13) yields 

( oo I 

{j+1) fH = r J (w - w b y i + E E U) 9m( 

{ 1=1 k=0 

j = 0,1,2,..., 

with the coefficients ^gn related to the coefficients ^gki by the formula 



(B.16) 



oo I 



1=1 k=0 

oo I 



oo I 



2a 



1=1 k=0 



1/2 



(B.17) 



1=1 k=0 



To get this formula, (B.16) with (B.17), we used again (B.8b) and we also set x = (w—w b )' y , 
y = (w — Wb) 2 ^^ ■ A closed form expression of the coefficients ^gui in terms of the 
coefficients ^)g kl is not available, but there is no difficulty in principle to compute them 
sequentially, by expanding the left-hand side of (B.17) in powers of x and y. 
It is then clear that (B.16) can be integrated to yield (via (B.6) and (B.7b)) 



(i+i) 



1 + E E (™ - ^) fc7 (™ - ^) 2i(1 ~ 7) 



1=1 k=0 



with 



ij+1) 9ki = (1 + k + 21a)- 1 W 



9kl- 



(B.18) 



(B.19) 



The consistency of (B.18) with (B.13) validates the ansatz (B.13) for all (nonnegative 
integer) values of the index j, and by assuming that this continues to hold in the limit 
as this index diverges, j — > oo, and that in this limit (w) yields the solution / (w) of 
(B.3a), as indicated by a comparison of (B.3a) with (B.5), we infer for / (w) the expression 



C oo I \ 

f(w) = 7 -7( w _ Wb y l + Y^Miw - w b p(w - w b ) 21 ^ , 
I 1=1 k=0 J 



3= 0,1,2,. 



(B.20) 



Here the coefficients gki are of course defined in terms of the above iteration, g^\ = ^°°^gkh 
and they depend only on the "coupling constant" a (and of course, as indicated by the 
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notation, on the two indices k and I), and we of course assume that, for w sufficiently 
close to w b , namely for \w — w b \ sufficiently small, the infinite sum over the index / in the 
right-hand side of this expression, (B.20), does converge. 

Of course we could have directly started from the ansatz (B.20) and then verified its 
consistency with the nonlinear ODE (B.3a) - but the iterative approach appeared to us 
preferable in order to explain here how one arrives at such an ansatz. 

Let us now proceed and discuss, always in the context of mechanism (i), the comple- 
mentary case to (B.7a), namely let us now assume that 

Re(a)<0, (B.21a) 

entailing of course 

Re(/3)>0. (B.21b) 
It is now convenient to rewrite the iteration formula (B.5) in the form 

y +1 )/'H=±|l+[ (j ' ) /H] } , 3 = 0,1,2,..., (B.22) 

and to solve it - consistently with (B.4a), again with the initial conditions (B.6)- by 
starting from the zeroth-order assignment 

^f(w)=±(w-w b ). (B.23) 

We then get 

«/' (w)=±{l + [±(w- w b )f} 1/2 , (B.24a) 

h = ± ji + £ ( 11 i ) [± (™ - ^} > ( B - 24b ) 

entailing, via (B.6), 

Wf (w) = ±(w-w b )^l + f2 (1 + Wr 1 ( 1 { 2 Y±(w- w b )f J . (B.25) 

We now set, as suggested by this formula, (B.25), 

®f (u,) = ±(«,-^){l + f {3) 9l [± (u> - Of J , j = 0, 1, 2, . . . , (B.26) 

and verify the consistency of this ansatz, (B.26), with the iteration (B.22). Indeed, by 
proceeding as above, we get 

H = ± (w - w b ) jl + k' +1 > 5i [± («, - «; b )f J , (B.27) 
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with 

V + Vg l = (l + l/3)- 1 Wgi (B.28) 



and 



l + x 



i + £ 



1=1 



1/2 

' = 1 + £°'W- (B.29) 



So, in analogy to the preceding case, we infer now the following expression for the solution 
of (B.3a) with (B.4a) and (B.21): 

/M = ±( w -«*){i + |> [± (™ - «>&)] ,/9 1 • ( B - 3 °) 

The significance of the choice among one or the other alternative for the ± sign is explained 
in Section 4. 

Let us now turn to the consideration of mechanism (ii), see (B.4b). In this case the 
assumption (B.7a) is inconsistent with (B.4b), because (B.4b) with (B.7a) entails, via 
(B.3a), that f'(wb) = ±1 , while a finite value of f'(wb) is clearly inconsistent with 
(B.4b). It is moreover clear that the assumption (B.7a) is consistent with (B.21) only if 
this inequality is strengthened to read 

Re(a)<-1, (B.31a) 

entailing (see (B.8a)) 

Re (7) < 0. (B.31b) 

Then the appropriate representation of the solution is given by (B.20), which is then 
clearly consistent with (B.4b) due to (B.31b) (and note that all the terms in the sums in 
the right-hand side of this representation, (B.20), become vanishingly small for w « w^, 
again thanks to (B.31b), because the index k is limited not to exceed the index I). 

This concludes our analysis, which validates the results reported in Section 4. 

Let us end this Appendix with a final remark, based on (B.l), which we conveniently 
rewrite in the following form (see (4.4a)): 

(C') 2 -BC 2a = V 2 , (B.32a) 

where of course V is given in terms of the initial data by (B.lb) and B by the formula 
(see (B.lb) and (4.1b)) 

B = ~< (0) C2 (0) [Ci (0) - C2 (0)] 2a . (B.32b) 

Assume now that the coupling constant a is a negative integer, a = —p, so that (B.32) 
reads 

(O 2 + < (o) a (o) [Ci (o) - c 2 (or 2p c 2p = y\ (B.33) 
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and assume moreover that the initial data, Ci (0), C2 (0), ([ (0), Q' 2 (0), are all real, so that 
this equation of motion, (B.33), and its solution, are also real (for real r), and the solution 
of the equations of motion (1.7) is as well real for real r. It is then clear from (B.33) (which 
can then be interpreted as the energy conservation formula for a one-dimensional particle 
in a power-law potential which vanishes at the origin and diverges at large distance) that, 
if the two initial "velocities" ([ (0), (' 2 (0) have the same sign, so that the potential energy 
(i.e., the second term in the left-hand side of (B.33)) diverges at large distances to positive 
infinity, then the solution ( (r) of this equation (B.33) is, for real r, always periodic, albeit 
with a period which depends on the initial data. In this case, however, the motion of the 
two "particles" Ci (t) and C2 (t) will not be periodic, due to the drift associated with the 
uniform motion of their center of mass (note that, if the two initial "velocities" (0), 
Q 2 (0) have the same sign, the center of mass will indeed move). Let us emphasize that 
this discussion refers to the solution of the equations of motion (1.7) in the special case 
in which they can themselves be interpreted to represent the (one-dimensional) motion of 
real particles. For explicit examples of this phenomenon see the treatments in Section 4 
of the two-body problems with pair coupling constants a = —1, a = —2 and a = —3, 
respectively. 



Appendix C 

In this Appendix we outline the treatment that justifies the main results of Section 5. 
Hence we investigate the behavior of the (generic) solution ( (r) of (1.7) in the neighbor- 
hood of its singularities. In analogy to what we did in Appendix B our basic approach is 
to replace (1.7) with the following recursive set of equations 

N 

^C = 2 £ a nm ^C'n U) C/( U) Cn- U) Cm), J = 0, 1, 2, . . . , (C.l) 

m=l; mjtn 

and to solve this recursive sequence of equations by starting from an appropriate assign- 
ment for the zeroth-order term (t) so that it capture the leading behavior of the 
solution £ (r) in the neighborhood of the singularity whose nature we wish to ascertain. 
Let us emphasize that this procedure is only supposed to work in the neighborhood of 
the value, r = t\,, at which a singularity occurs, and its effectiveness is predicated upon 
the successful identification of the zeroth-order assignment (°) £ (r) , as demonstrated by 
the subsequent preservation of the leading part of its behavior for r « r\, throughout the 
iteration (see examples below). We restrict our treatment to a demonstration of this fact; 
to qualify as a complete proof our analysis should be complemented by two additional, 
synergistically related, results (which will instead be taken for granted here): (i) a proof 
that the resulting representation of ( (r) as an infinite series (see for instance (5.4), (5.6) 
and (5.9)) does converge in a sufficiently small neighborhood of the singular point r = tj,, 
and (ii) a proof that the sequence produced by the iteration (C.l) (supplemented by ap- 
propriate "initial conditions" at every iteration cycle: see examples below) does converge 
to a solution ( (t) of (1.7), 

lim [^Cl = C (r) . (C.2) 
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Throughout this Appendix C we use the notation of Section 5 without reporting anew 
the corresponding definitions. Our first task is to justify (5.4). To this end we set 

(0) Ci (r) = b + c (r - r b y + v(T-r b ), (C.3a) 
(0) C2 (t) = b c (t r b f + « (r - r 6 ) , (C.3b) 
(0) Cn (r) = 6 n + v n (r - r b ) , n = 3, . . . , TV. (C.3c) 

From this assignment and (C.l) (with j = and n = 1) we get 

WCf = 2a[c 7 (r - nY'- 1 + v}[- n (r - n)^ 1 + v]/[2c(t - nf] 

N 

+ 2 ai m [c 7 (r - Tfe) 7 " 1 + v]v m /[b - b m + c(r - r fe ) 7 + (« - u to )(t - r 6 )] (C.4a) 

m=3 

= _ a7 2 c(r _ ^7-2 (1 _ v 2 (c7) -2 (r _ T6) 2-2 7] 

+ 2c 7 (r - ny-^i + Kc7) _1 (r - n) 1 ^} 

N 

x £ a lm z; m (6 - b^- 1 {1 + (6 - M^Kr - r b y + (v- v m )(r - n)}}' 1 . (C.4b) 

m=3 

The motivation to re-write (C.4a) in the form (C.4b) is of course based on (5.2b), which 
is as well essential for the subsequent developments. 
We now use the following formula: 



fc=0 z=o 

where 

' k 



1 , -k\/[l\(k-l)\] (C.5b) 

is the binomial coefficient, which is hereafter understood to vanish if I or (A; — Z) is a 
negative integer. This formula, (C.5), is of course only applicable provided 

|e + <5|<l. (C.6a) 

Via (C.5) we rewrite (C.4b) as follows: 

(1) Ci = 7(7 " l)c(r - T h y- 2 + (av 2 /c)(r - r b )^ 

oo k 

+ [l + v (C7)' 1 (r - n)^} £ £ /£> (r - t^" 1 * 1 ^- 1 , (C.4c) 



where 

N 



l J ^2 a ^rnV m {b - bm)' 1 
' m=3 



x [-c/(b - b m )] k - l [-(v - v m )/(b - b m )} 1 , (CM) 
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and we also took advantage of the relation 

7 - 1 = -07 (C.7) 

entailed by the definition (5.2a) of 7. Note that the step from (C.4c) to (C.4d) is justified 
provided the following condition holds (which corresponds to (C.6a)): 

\c(r - r b y + (v- v m )(r - r b )\ < \b - b m \ , (C.8) 

which is certainly true (see (5.2b)) provided r is sufficiently close to r b , namely |r — T b \ is 
sufficiently small. 

We now integrate twice (C.4c), adjusting the two integration constants so that, for 
r ~ r b , the three leading terms of (t) coincide with (r), see (C.3a) (note that 
the first correction term next to the constant one, namely the second one of the three, is 
automatically OK) and we thereby get 



(i) 



Ci (r) = 6 + c (r - r b y + V {T-T b ) + av 2 [c (2 - 7) (1 - 7)]- 1 (r - n) 2 ^ 



00 00 



+ EE (r ] $ (r ~ nY (r - n) k \ (C.9a) 
j=i k=i 

with 

= u + kjy 1 (j - 1 + ^r 1 /iiUj-i 

+ v (C7)" 1 \j + l + (k- 1) 7]" 1 [j + (k - 1) 7P 1 f$ k -2,y (C9b) 
In a completely analogous manner we get 

( 2 )Ci (r) = 6 - c (r - + v (r - r 6 ) - at; 2 [c (2 - 7) (1 - 7)]^ (r - r b ) 2 ^ 

00 00 

+EE ( r - ( r - T ^ > ( c - 10a ) 
i=i fc=i 

with 

( \ ( ? = (j + fry) -1 (j - 1 + fry)" 1 

- v (C7)- 1 [j + 1 + (k - 1) j}' 1 [J + (k~l) I}' 1 f™ k _ 2d , (C.lOb) 



where 



fid = ~ 2cr y ( / ) E a 2mVm(b - b m ) 1 



m=3 

X [c/(6 - 6m)] fc "'[-(^ " " bm)] 1 . (C.lOc) 

Also analogous is the derivation of ( 1 )^ n (r), n = 3, . . . , TV: 
(1) Cn(r) = 6„ + f„(r-T 6 ) 

OO OO 

+ E E ( T " ^ ( T - r ^) fc7 > n = 3,...,N, (C.lla) 

j=l fc=5ji 
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= l(j + in) (3-1 + fry)]" 1 1 (i - s n ) K/ (b n - b)] [ C / (b n - b)] k 

x [-(„„ - z,)/(6 n - 6)P'- 2 [^ + fc(j - iy l l(v - v n )\ 
x(^J- 2 ) 2 ^a,(-l)(^ 

^ ' s=l,2 

^ -1 -1 ] 

+ 6 kfi 2^2a nm v n v m (b n -b m ) 1 [- (v n - v m ) / (b n - b m )} 3 1 \. (C.llb) 

m=3 ) 

At this point we introduce the ansatz 

(r) = b+ (-I)*- 1 c (r - r b y + v(t- n) 

OO oo 

+ E E ( ^l(--^ + ' 7+m(1 - T) ' - = 1,2, (C12a) 

fc=l Z,m=0;Z+m>l 

(i) Cn (r) = b n + v„ (r - n) 

OO OO CO 

+ E E E ^ffii - r b r^ m ^ , n = 3, . . . , N, (C.12b) 

fc=l /=5 fe i m=0 

which is clearly consistent (for the appropriate choice of the coefficients g^ m ) with both 
(C.3) (j = 0) and (C.10,11) (j = 1), and we show that it is moreover consistent with 
the iteration (C.f), namely that insertion of this ansatz, (C.12), in the right-hand side 
of (C.l) yields, after appropriate expansions and integrations, its validity at the iteration 
order j + 1. To this end one must of course use the formal expansion formula 

( oo ^| — 1 oo 

|l + E^[ =1 + E^e fc , (C.13a) 
I k=i J k=i 

where we assume that the coefficients are given and the coefficients are obtained from 
them (clearly this can be done sequentially starting from the lowest powers, although a 
closed form expression is not generally available). Of course this formula can only hold 
under the restriction 

|e| < E, (C.13b) 

where the constant E is required to be small enough so that all the sums in (C.13a) 
converge (and in order that this be possible appropriate restrictions must hold to begin 
with on the behavior of the coefficients g k when their indices diverge, to guarantee that 
the sum in the left-hand side of this formula, (C.13a), do converge when (C.13b) holds). 

It is easy to convince oneself that all one in fact needs to prove is that the following 
expansions, 

v\ 9 (t) = b + (-i)(-v c (t - n y + V(T - T b ) 

+ 0(\(t- n )\ 1+ \ \(r - T b )\ 2 ^) , 8 = 1,2, (C.14a) 
U) Cn(r) = b n + v n (T-r b ) + o(\(T-r b )\ 1+ ^), n = 3,...,N, (C.14b) 
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which are clearly entailed by (C.12) (and are of course valid for r « rfc), are consistent 
with the iteration (C.l), namely that their validity at the order j, see (C.14), entails their 
validity at the order j + 1. Indeed this result (once proven) entails that the iteration 
process keeps producing only the powers of (r — T&) that appear in the right-hand side 
of (C.12) (and no logarithms!), with the coefficients g^ m uniquely determined via 

(C.l) by the coefficients ^g^ m and the other constants appearing in the right-hand side 
of (C.12). 

To prove the consistency of (C.14) with (C.l) we now note that these formulas, (C.14), 
entail 



C*)£(r) = 7 (-l)(- 1 ) c ( r - n)^ 1 + v 

+ 0(\(T-r b W\,\(T-T b )\ 1 -'), a = 1,2, 
W£(T)=v n + 0(\(r-n)\' r ), n = 3,...,N, 



(C.15a) 
(C.15b) 



hence 



and 



{j) C'i (r) W CUr) = -7 2 c 2 (r-r^- 1 )+, 2 

+ o(\(T-n)\ 2rr - 1 A(r-n)\), (C16a) 

{j) C(r)^C n (r) = (-l)(-D 7CUn (r - rtf' 1 +vv n + (|(r - r 6 )n(r - r^) , 

s = l,2, n = 3,...,N, (C.16b) 
(J) C (r) (j) C 00 = t>n«m + O (|(r - r b )D , n, m = 3, . . . , N, (C.16c) 



(j) Ci (r) - ^C2 (r) = 2c(r - r b f + O (|(r - r b )| 1+ ^ , |(r - r b )| 2 ^) , (C.17a) 
W C s (r)-Wc n (r) = 6-6 n + 0(|(T-r 6 )n, s = 1,2, n = 3,...,N, (C.17b) 
(i) Cn(r)-^C m (r) = 6 n -6 m + 0(|(T-r fe )| 1+7 ), n, m = 3, . . . , iV, (C.17c) 



hence 



U) Ci(r)- U) C2(r) 



(2c)-\T- n )-i 

x "l + O^T-r^KKr-r;,)! 2 ^) 



a) C s (r) - ^Ur)\ =(b- b n y l [1 + O (|(r - r b )P)] , 

s = l,2; n = 3, ...,iV, 
a) Cn(r) - (j) C m (r)] _1 = (6 n - b m )~ l [l + O (|(r - r b )| 1+ ^) 



(C.18a) 



(C.18b) 



n, m = 3, . . . , N. 



(C.18c) 
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It is therefore clear that (C.l) entails (via (C.16a,b) and (C.17a,b)) 

&" +1 >C (r) = 2a [- 7 V (r - r b ) 2 ^ + v 2 + O (|(r - ti)^" 1 , |(r - r 6 )|) 
x (-I) 8 " 1 (2c)- 1 (r - r 6 H [l + O (\(t - r b )\ , |(r - ^^(i-y)' 



A? 



+ 2 a sm [(-l) (s " 1} 7™n (r - r,)^ 1 + VVn + O (|(r - r b )P , \(r - n)] 1 ^) 

m=3 

xib-bny^l + Oiftr-nW)], s = l,2, (C.19a) 

namely 

{i+1) Cs(r) = (-irS(7 - 1)(t - r 6 )^" 2 

+ 0(|(r-r 6 )r,|(r-r 6 r 1 |) 1 



s = 1,2. 



(C.19b) 



Integration (with the appropriate choice of integration constants) of this equation, (C.19b), 
yields 



^ +1) Ut) = b + (-l)^" 1 )^ - r b y + v(t - r b ) 

+ 0(\(r-T b )\ 1+ \\(T-n)\ 2 -^, 8 = 1,2, 

and the consistency of this expression, (C.20), with (C.14a) is plain. 
It is likewise clear that (C.l) entails (via (C.16a,c) and (C.17a,c)) 



(C.20) 



VVr, 



(r) = 2^ [a ns [ (-l)C- 1 ) 7 cv n (r - nf' 1 + ^ 

8=1,2 

+ 0(\(r- r b )\\ \(r - n)^, \(r - r b )\ 2 ^ 1 ) ] (b n - b)' 1 [1 + O (|(r - r b )P)] } 

N 

+ 2 £ {a nm M m + O (|(r - r b )| 7 )] (6 n - ^[l + O (|(r - r 6 )| 1+7 )]} , 

m=3 

n = 3,...,iV, (C.21a) 

namely 

(r) = O (|r - nf' 1 ) , n = 3,...,N, (C.21b) 

and the integration (with the appropriate choice of integration constants) of this equation, 
(C.21b), yields 



0+1) 



Cn (r) = 6 n + v n (r - r b ) + O (|(r - r 6 )| 1+7 ) , n = 3 



....,iV, 



(C.22) 



which is clearly consistent with (C.14b). (Note however that the estimate (C.21b) can be 
replaced by more stringent estimates if a n \ = a n 2 and v n ^ 0, or if v n = 0; indeed it is 
easily seen from (C.21a) that the last term in the right hand sides of (C.14b) and (C.22) can 

be replaced by O ( |r — T b \ 2 ) if a n \ = a n 2 and v n ^ 0, and by O ( |r — r b | 2+7 , |r — Tbj 3-7 ) 



if v n = 0. This entails in subsequent equations some corresponding changes that the 
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reader will easily ascertain, but which do not modify the main conclusions, see Section 5, 
as regards the nature of the branch point of the solution ( (r)). 

This concludes our treatment meant to justify the validity of the representation (5.4). 

Let us now proceed to justify the ansatz (5.6), namely the compatibility of the expan- 



sion 



C 1 (r) = b + c(r-r b )^ + £ $ (r - r b f^ , 

k,l=l;k+l>3 

oo oo 

C2{T) = h+V2{T - Tb) - c{T - Tb f+p + Y^ £hr-nt +ip . 

k=l l=28 kl 



,N, 



(C.23a) 
(C.23b) 
(C.23c) 



Cn (r) = b n + v n {T-r b ) + Y J Yl 9m (t - nf^ , n = 3, . 

k=2 1=0 

where 

(3 = -2a, Re(/?)>0, (C.23d) 

with the evolution equation (1.7) satisfied by the solution C( r )- Of course this ansatz, 
(C.23), is again only applicable for r « r b , namely in the neighborhood of the value r b of 
the independent variable where the solution ( (r) features a "two-body collision" , charac- 
terized, see (C.23), by the equality and inequalities 

Ci (n) = (2 (n) + Cn (n) , n = 3, . . . ,N. (c.24) 

And clearly in that neighborhood this ansatz, (C.23), entails 



Ci (r) = b + c (r - T h ) 1+P + O (Jr - T b \' +f " , |r - r b 
C2 (r) = b + v 2 (r -n)-c(T- T b ) 1+(B + O (|r - r 6 | 2 , |r - r 6 | 1+2/3 ) 
Cn (t) = b n + w n (r - T b ) + O - r b | 2 ^ , n = 3, . . . , 



|2+/3 



1 1+2/3 



AT. 



(C.25a) 
(C.25b) 
(C.25c) 



Here we limit our treatment to show that this behavior is compatible with the evolution 
equation (1.7), which is the key point of our argument (see above - it is instead left for the 
diligent reader to repeat, in close analogy to what we did before, the iteration argument 
that leads to the ansatz (C.23)). 

Clearly the expansions (C.25) entail 



C' 1 (r) = (l + P)c(r-T b f 



l + O Ir-rJ 



C2 (t) = V 2 



C' n (r) = v n [l + 0(\r-T b \ 



1 + 0[\t- n\ , \t - r b \ 
,\T-T b f)\ 



hence 



n = 3,...,N, 
1 + o(\t- 



C[(r)C(r) = (l + (3)cv n (T-r b ) 13 

n = 2,...,N, 
(2 (r) C (r) = V2V n [l + O (|r - r b \ , |r - ^l") 



n\ 



-nf) 



C (t) C' m (r) = v n v m [1 + O (|r - T b \ 



n, m 



n = 3,...,N, 
...,N, 



(C.26a) 

(C.26b) 
(C.26c) 



(C.27a) 
(C.27b) 
(C.27c) 
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as well as 

Ci (r) - C2 (r) = -v 2 (t - r b ) [l + O (|r - r 5 | 1+/3 )~ 

C.(r)-Cn(r) = (6-6„)[l + 0(|r-7i|)], a = l,2, n = 3,...,iV, 
Cn (r) - Cm (r) = (b n - b m ) [1 + O (|r - T 6 |)] , n, m = 3, . . . , N, 



hence 



[Ci (r) - C2 (r)]- 1 = -wj 1 (r - n)' 1 [l + O (|r - r 6 | 1+ ^ 

[Cj (r) - Cn (r)]- 1 = (b - b n )' 1 [1 + (|r - r b \)\ , 

s = l,2, n = 3,...,N, 
[Cn (r) - C™ (r)]" 1 = (6„ - [1 + (|r - r 6 |)] , n, m = 3, . . . , TV. 



(C.28a) 

(C.28b) 
(C.28c) 



(C.29a) 

(C.29b) 
(C.29c) 



It is now easy, using these formulas, to verify the consistency of (C.25) with (1.7). 
Indeed for n = 1 (1.7) yields (via (C.27a) and (C.29a,b)) 



Cf = -2ac(l+0) (T-r b ) 



-1+13 



X 1 + O (|r - 7i| , \t - Til")] + O (|t - nf) 
namely, via (C.23d), 

C'{ = c0 (1 + /?) (r - r b r 1+/3 + O (|r - rfcl' 3 , |r - r b |" 



/ 3 



(C.30a) 



(C.30b) 



which can be immediately integrated (with an appropriate choice of integration constants) 
to yield (C.25a). Likewise for n = 2 (1.7) yields (via (C.27a) and (C.29a,b)) 



l + 0(\T-T b \,\r-r b f)] +0(1), 



C£ = 20C(1 + /J) (T-T b )- 1+/3 

namely, via (C.23d), 

= -c(3 (1 + 0) (r - r,)" 1 ^ + O (|r - r,!" 1 ^ , l) , 



(C.31a) 



(C.31b) 



which can be immediately integrated (with an appropriate choice of integration constants) 
to yield (C.25b). And finally for n = 3, . . . , N (1.7) yields (via (C.27b,c) and (C.29b,c)) 



C = 0(1), 



(C.32) 



which can be immediately integrated (with an appropriate choice of integration constants) 
to yield (C.25c). 



Appendix D 



In this Appendix D we introduce the monic polynomials 11^ (C, t) respectively Pjv (z, t), of 
degree TV" in the variable C respectively z, the zeros of which coincide with the coordinates, 
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Cn = Cn (t) respectively z n = z n (t), of the particles moving according to the A-body 
problems (1.7) respectively (1.5), 



TV TV 



% (C, r) = C" + £ 7m (r) C^"™ = II [C - Cn (r)] , (D.l) 

m=l n=l 
TV TV 

P N (z, t) = z N +J2 c m (t) z N ~ m = l[[z-z n (t)} . (D.2) 
m=l n=l 

Via such a position a one-to-one correspondence is introduced among the (ordered) sets 
of coefficients (r) , m = 1, . . . , A} respectively {c m (t) , m = 1, . . . , A} and the (un- 
ordered) sets of zeros {( n (r) , n = 1, . . . , A} respectively {z n (t) , n = 1, . . . , A}. This 
correspondence turns out to be extremely convenient [5, 2] in the integrable case in which 
all the coupling constants a nm in (1.5) and (1.7) equal unity, 

a nm = 1. (D.3) 

Indeed in this integrable case these polynomials satisfy extremely simple linear PDEs, 
namely (note that throughout this Appendix D subscripted variables denote partial dif- 
ferentiations) 

njv, rT = 0, (D.4a) 
Pn,u = iuPN,t, (D.5a) 

entailing of course 

7™ = 0, 7m (r)= 7m (0)+7^(0)r, (D.4b) 
c m = iujc m , c m (t) = c m (0) + c m (0) [exp (iut) - 1] / (iu) . (D.5b) 

It is this connection that lies behind the resolvent formula (1.17). 

In this Appendix D we investigate the extent to which an analogous approach is useful 
to treat the A^-body problems (1.7) and (1.5) in the (generally nonintegrable) case when 
the coupling constants a nm are still all equal, but they need not be unity, 

a nm = a, (D.6) 

with a an arbitrary (possibly complex) constant. 

Hereafter we focus on the A-body problem (1.7), and only at the end, by taking 
advantage of our by now usual trick (see (1-6)), we mention the analogous results related 
to the A-body problem (1.5). We also drop, whenever this is notationally convenient, the 
subscript N attached to the symbols II and P (namely hereafter il/v = II, P/y = P). 

Our results hinge on the following basic result, the proof of which is relegated to the 
very end of this Appendix. 

Lemma D.l. If the A zeros Cn ( T ) evolve according to (1.7) with (D.6), the (monic) 
polynomial II ((,r) of degree A in Q, see (D.l), satisfies the following nonlinear PDE: 

u TT [n c ] 2 + (l - a) |n cc (n r ) 2 - n c [(n T ) 2 ] J = *n, (d.7) 
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where VP = VP (£, r) = ^2N-4 (C r ) *s a polynomial of degree 2N — 4 in the variable (, 

2N-A 

*(C,r)= J] ^MC" - *- (D.8) 

fe=0 

Note that this PDE, (D.7), does not explicitly feature the parameter N. 
Let us note that (D.l) entails the identification 

N 
n=l 

so that the coefficient 71 (r) clearly, see (1.7), satisfies the equation 

7" (r) = 0. (D.lOa) 

Hence U TT (£, r) is a polynomial in ( of degree N — 2 (see (D.l)). It is therefore clear, see 
(D.l), that the two sides of the PDE (D.7) are polynomials of degree 3N — 4 in £, hence 
this PDE entails 3N — 3 equations to be satisfied by the N coefficients -y n (r), see (D.l), 
and the 2N — 3 coefficients ipk (t), see (D.8). Hence the number of equations matches the 
number of unknowns (3N — 3 = N + (2N — 3)). Moreover these equations are clearly all 
linear and purely algebraic (not differential!) for the 2N — 3 coefficients ip^ (r), see (D.7) 
and (D.8). It is therefore in principle possible to eliminate the 2N — 3 unknowns ipk (r) 
and to thereby obtain a system of only N, generally nonlinear, second-order ODEs for the 
N unknowns 7„ (r), or rather a system of N — 1, generally nonlinear, second-order ODEs 
for the N — 1 unknowns 7„ (r), n = 2, ... ,7V, since 71 (r) shall turn out to satisfy the 
linear second-order ODE (D.lOa), which of course entails 

7i(t)=7i(0)+7i(0)t. (D.lOb) 

A systematic algorithmic procedure ("Euclid's method") to eliminate the coefficients 
of the polynomial ^ from (D.7) and to thereby obtain the (generally nonlinear) second- 
order ODEs satisfied by the coefficients 7„ (r) can be based on the observation (which is 
actually instrumental, see below, to obtain the PDE (D.7)) that the right-hand side of 
(D.7) vanishes for ( = ( n , see (D.l). One can therefore lower the degree of the left-hand 
side of (D.7) by applying sequentially the substitution 

N 

C N = ~J2 7™( N - m (D.H) 
m=l 

to the term of highest degree, or to all terms of degree higher than N. After applying this 
procedure 2N — 3 times the left-hand side of (D.7) becomes a polynomial in £ of degree 
N — 1, since each application lowers by one unit the degree of the polynomial. But since 
the polynomial of degree N — 1 thus obtained must vanish at the TV points ( = Cn (see 
(D.7)), it vanishes identically, namely all its N coefficients vanish: and this requirement 
yields precisely the N equations satisfied by the iV coefficients 7„ (r) (hence one of these N 
equations will be precisely (D.lOa)). For instance for N = 2 one thereby obtains 

ii = 2 (I -a) (4 72 - 7?)" 1 (7I - TWiTi + 1%) , (D.12) 
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of course with 71 (r) given by (D.lOb), and for N = 3 one obtains 

H = 2(1 - a)5~ l [ 7 f (-67! 7l 2 - 47x7273 + 7! + 9 7 I) + I2 (-3ti73 + 7 2 2 ) 

+ 7 3 2 (7 2 - 372) + 7i7 2 (-7172 + 4 7?73 - 37273) + 27W3 (-37173 + ll) 

+ 7 2 73 (-37371 + 7 2 2 )], (D.13a) 

7 3 ' = 2(1 - a)^ 1 [7^3 (-37173 + 7 2 2 ) + 7 2 2 73 (7? " 37 2 ) + l'i (7? " 4 7i72 + 9 73 ) 
+ 7i7273 (-7172 + 973) + 27(7373 (7? - 372) 

+ 7 2 7 3 (-37173 - 7?72 + 4 72 2 ) ] , (D.13b) 
6 = 27 7s 2 - I8737271 + 4737? + 47! - 7 2 2 7i 2 , (D.13c) 

of course with 71 (r) given again by (D.lOb). 

Note the consistency of these equations with the general result 7" = 0, valid for arbi- 
trary N in the integrable (a = 1) case [5, 2]. 

The diligent reader will also verify the consistency of these ODEs, (D.12), (D.13) and 
(D.7), with the solutions in the free (a = 0) case, when of course, for all n = 1, . . . , N with 
arbitrary N, ( n (r) = ( n (0) + Cn (0) T ( see an d (D.6), with a = 0), and the coefficients 
7 n (r) have the corresponding expressions entailed by the well-known explicit formulas 
expressing the coefficients of a polynomial in terms of its zeros, see (D.l). The diligent 
reader will also obtain the solutions of (D.12) that correspond to the two-body problems 
treated in Section 4, and the solutions of (D.13) respectively (D.7) that correspond to the 
similarity solutions (with N = 3 respectively arbitrary N) of Section 3. 

The results of [11] entail the conjecture that, if a = —1/2, the solutions of the two non- 
linear coupled nonautonomous ODEs (D.13) with (D.lOb) possess the ("super Painleve") 
property to only feature solutions which are entire functions of the independent variable r. 

Let us now report tersely the analogous results relevant to the equations of motion (1.5) 
rather than (1.7). To this end we set 

U n (C,t) =P N (z,t), ( = z, t = [exp (iwt) - 1] / (iw) , (D.14) 

as suggested by (1.6), and we then note that the nonlinear PDE (D.7) can now be re- 
written as follows: 

(P tt - iujPt) [P z ] 2 + (1 - a) [P zz {P t f - P z [(P t ) 2 ] } = QP- (D.15) 
Here we introduced via 

Q = Q(z,t) = exp (-2iwt) * (C, r) (D.16) 
the polynomial Q (z,t) of degree 2N — 4 in the variable z, 

2N-A 

Q(z,t)= Qk{t)z N - k . (D.17) 

fc=0 

We do not elaborate any further on the polynomial solutions of the PDE (D.15), since 
with trivial modifications what is written above after (D.7) remains applicable, and it 
is as well quite easy to apply the change of (independent) variable j n (r) = c n (i), see 
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(D.14) and (D.1,2) in order to obtain from (D.lOa), from (D.12) and from (D.13) the 
corresponding ODEs for the coefficients c n (t). We only note that we expect (D.15) to 
possess generally a lot of solutions that are completely periodic with period T, see (1.2), in 
the (real) independent variable t. And in particular it is natural (see above) to conjecture 
that, if a = —1/2, all polynomial solutions P (z, t) of degree N = 3 of the nonlinear PDE 
(D.15) are completely periodic with period T in the (real) independent variable t. 

Finally, let us obtain (D.7) from (D.l) and (1.7) with (D.6). 

Differentiation of (D.l) with respect to r yields 

N N N 

J2^'m(r)C N - m = -J2C(r) II [C-O(r)], (D.18a) 

m=l n=l 1=1; tyn 

and a further differentiation yields 

N N N 

E^K iV " m = -ECn(r) II [C-OW] 

m=l n=l 1=1; l^n 

N N N 

+ E C/(r) II [C-Cfc(r)]. (D.18b) 

n=l l=l;lyin k=l;k^l,n 

We now set C = Cn 111 (D.18a) and we thereby get 

N 

C (r) = ~ E -L (r) [Cn (r)] N ~ m /U ( [( n (r) , r] , (D.19) 

m=l 

where we used the relation 

N 

^[Cn(r),T}^dU(C,r)/dC\ (=UT) = J] [Cn(r)-Ci(r)]. (D.20) 

1=1; l^n 

Likewise, by setting ( = ( n in (D.18b) we get 

N 

-C(t) + 2C(t) E C/(r)/[Cn(r)-0(r)] 

l=l;l^n 

N 

= E ^{r)[Ur)] N - m /n C [Cn(r),r], (D.21) 

m=l 

where we used again (D.20). 

We now use the equations of motion (1.7) and we thereby get 

N 

K(r) E (l-a nl )Ci(r)/[Ur)-Cl(r)] 

1=1; tyn 
N 

= E l'M[Ur)} N - m /U c [Cn(r),r], (D.22) 

m=l 
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and via (D.19a) this yields 

N N N 

2 £ I'm (r) [Cn (r)] N - m E (i-MEiMldM]' 

m=l l=l;l^n k=l 

N 

x {[ Cn (r) - (r)] n c [0 (r) , r}}' 1 = £ 7m (r) [C„ (t)]""™ . (D.23a) 

m=l 

From now on we restrict consideration to the equal-coupling-constants case, see (D.6), 
and we re- write (D.23a) as follows: 

N N 

2 E ^m(r)lk(r)[Cn(r)] N - m £ (1 - a) [Q (r)]""* 

m,k=l l=l\l^n 

N 

x {[Cn (r) - (r)] n c [0 (r)^}}- 1 = E 7 m (r) [Cn (r)]^" 1 . (D.23b) 

m=l 

We now note that the rational function of (, 

F nl (C) = C N - k l(Cn-C)n(CT\ (D.24) 

where 11(C) is the polynomial of degree TV that has the TV zeros Cn (see (D.l)) and k is 
an arbitrary positive integer not larger than N, k = 1, . . . , TV, vanishes faster than C _1 as 
C — > oo, so that the sum of all its residues vanishes. Since this rational function has TV — 1 
simple poles (at C = Cm, m = 1, . . . , N , m ^ n) and a double pole (at C = Cn), this implies 
the identity 

N 

E c/^KCn-owor 1 

1=1; tyn 

= (N- AOC - *- 1 ^)]- 1 - (i/2)c^- fc n cc (c„)[n c (c„)]- 2 . (D.25) 

Via this identity (D.23b) yields 

N 

(1-a) E Vm ( T ) ^ ( T ) { 2 (TV - fc) [Cn (r)] 2iV — n c [C n (r) , r] (D.26) 

m,k=l 

AT 

- [Cn (r)] 2N - m - k U a [C n (r) , r] } = (H f (r) , r]) 2 E 7^ (r) [Cn (t)]""™ . 

m=l 

We now note that the definition (D.l) of the polynomial II (C, t) implies the relations 

N 

E (r) 7^ (r) [Cn (r)] 2JV — * = [U T (C, r)] 2 | c=Cn(r) , (D.27a) 

m,fc=l 

AT 

E 7 m (r) 7^ (r) 2(N — k) [Cn (r)] 2 ^^" 1 (D.27b) 
l 

= E 7m ( T )7fc ( r ) ( 2iV -m-k) [Cn ( T )] 2Ar - m - fe - 1 = |[n T (C,r)] 2 | I 



m,k=l 

N 
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where the first step in (D.27b) is of course justified by the possibility to exchange the 
dummy indices m and k. Hence from (D.26) and (D.l) we infer that the polynomial 



of degree 3N — 4 in £, vanishes at the N points £ = (, n , namely that it has the same zeros 
as the polynomial of N degree II, see (D.l), and this clearly implies the validity of (D.7). 
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